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ABSTRACT 

The  results  of  an  experiment  to  test  the  utility  of  satellite-observed 
data  for  enhancing  the  accuracy  and  detail  of  numerical  analyses  and  its 
effect  on  the  attendant  primitive  equation  (PE)  prognoses  is  presented. 
The  test,  carried  out  in  the  framework  of  Fleet  Numerical  Weather  Central's 
system,  employed  SIRS  data  and  various  quasi  objective  techniques  to 
establish  bogus  pressure/height  and  wind  data  at  sea  level/500  mb  over 
the  North  Pacific  Ocean  area.    The  resulting  satellite-modified  analyses 
for  five  synoptic  times,  in  the  period  10-12  March  1971, were  used  to  gen- 
erate PE  prognoses  from  1200  GMT  10  March  and  0000  GMT  11  March  1971. 
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tional  data  indicated  relatively  more  improvement  at  sea  level  than  at  500 
mb.    The  complex  results  are  interpreted  in  terms  of  the  initialization  of 
the  analyzed  fields,  known  behavioral  characteristics  of  the  PE  model, 
the  mechanics  of  the  modification  procedures  and  the  experience  level  of 
the  analysts  . 
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I  .      INTRODUCTION 

This  thesis  reports  on  a  preliminary  scheme  for  integrating  four 
methods  of  converting  satellite-observed  data  into  conventional  parame- 
ters for  incorporation  into  the  data  base  of  an  operational  numerical 
analysis  program.    The  four  methods  utilized  are:     1)  subjective  changes 
to  the  current  analyses  based  on  the  interpretation  of  video  and  infrared 
(IR)  weather  satellite  observations;    2)  cloud-motion  vector  winds  from 
satellite  time-lapse  photographs;   3)  statistically-derived  500  mb  heights 
from  satellite-observed  cloud  patterns;  and  4)  statistically  derived 
pressure-temperature  soundings  from  satellite  infrared  spectrometer  (SIRS) 

i"-\/^i  -\v^/»a       j~\  »-\  rr*  /"\vr  r  — i  +■  i  ^\tr\  <-! 

The  data  derived  from  these  methods  were  used  in  the  objective  re- 
analyses  of  the  Fleet  Numerical  Weather  Central,  Monterey  (FNWC), 
sea-level  pressure  and  500  mb  height  fields  at  five  synoptic  times  for  the 
period  10-12  March  1971.    Two  dates  were  selected  for  making  parallel 
numerical  prognoses,  with  and  without  satellite  data,  by  using  the  FNWC 
primitive  equation  (PE)  model. 


In  the  foregoing  discussion,  the  following  terms  are  used: 
UNMOD  ANAL  —  analysis  derived  from  conventional  data  only. 
MOD  ANAL       — analysis  derived  from  conventional  data  and  satellite 

observations . 
UNMOD  PROG-- prognosis  generated  from  the  UNMOD  ANAL. 
MOD  PROG       — prognosis  generated  from  the  MOD  ANAL. 
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Though  the  generation  of  quantitative  data  from  the  satellite  obser- 
vations was  not  accomplished  in  real  time,  the  scheme  followed  may  be 
adapted  easily  to  the  operational  schedule  at  a  major  computer  center 
which  has  access  to  satellite  observations. 
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II.      BACKGROUND 

The  major  meteorological  centers  rely  almost  entirely  on  conven- 
tional synoptic  data  (i.e.,  surface  observations,  rawinsondes,  aireps, 
etc.)  for  their  operational  numerical  analysis  and  prediction  schemes. 
The  question  of  incorporating  largely  qualitative  satellite  video  and  radi- 
ance observations  into  this  automated  computer  system  is  frequently 
referred  to  as  the  SINAP   (satellite  input  to  numerical  analysis   and  pre- 
diction)  problem  [8]    . 

Conventional  sources  of  data  are  frequently  scarce  over  the  Pacific 
and  Atlantic  Oceans .    Here,  satellite-observed  cloud  patterns   attain 
their  best  development,  and  radiance  measurements  have  theii    leatl 
contamination  owing  to  the  absence  of  hot  or  high  terrain,    snow  or  ice 
reflection,  and  so  forth.    Consequently,  the  thrust  of  satellite  meteorology 
this  past  decade  has  been  to  quantify  satellite  information  from  these 
oceanic  regions,  particularly  the  North  Pacific  Ocean,  because  the  poor 
conventional  data  base  there  can  lead   to  erroneous   numerical  predic- 
tions in  terms  of  accurate  weather  forecasts  for  North  America. 

Four  methods  for  attacking  the  SINAP  problem  have  received  strong 
emphasis.    The  general  background  of  these  methods  follows. 
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A.  SUBJECTIVE  ANALYSIS  CHANGES  FROM  THE  INTERPRETATION  OF 
VIDEO  AND  INFRARED  WEATHER  SATELLITE  OBSERVATIONS 
Subjective  changes  to  either  manual  or  objective  analysis  have  been 

undertaken  since  the  first  meteorological  satellite  video  pictures  were 
received  in  1960.    Global  video  coverage  was  realized  with  the  launching 
of  TIROS  LX  in  1965.    Since  1964,  infrared  sensors  have  provided  displays 
of  cloud  patterns  similar  to  the  video  pictures ,  but  with  significant  diag- 
nostic differences .    Global  infrared  nighttime  coverage  was  made  opera- 
tionally available  with  the  launching  of  ITOS   I  in  1970. 

In  general,  the  synoptic  meteorologist  interprets  the  position  or 
intensity  of  important  weather  features  from  a  visual  study  of  current 
satellite  pictures.    He  then  changes  the  operational  analyses  by  either 
manually  re-analyzing  the  weather  charts  or  by  fabricating  "bogus" 
pressure,  temperature,  or  wind  values  to  insert  into  a  computer  re-analysis 
scheme. 

Numerous  studies  report  the  substantial  progress  made  in  the  inter- 
pretation of  synoptic-scale  cloud  patterns  and  the  incorporation  of  the 
satellite  data  into  subjective  or  objective  (numerical)  analysis  schemes. 
The  reader  is  especially  directed  to  Ref s .   1,  7 ,   12,   14,   17,  and  26  for 
video  interpretations  and  to  Ref.   15  for  infrared  imagery  interpretations. 
Terminology  developed  in  these  references  will  be  used  in  this  thesis. 

B.  CLOUD-MOTION  VECTOR  WINDS 

Since  the  earliest  satellite  pictures  of  cloud  patterns,  meteorolo- 
gists have  attempted  to  obtain  wind  estimates  from  such  features   as 
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cirrus  "blow-off"  from  cumulonimbus  clouds,  cirrus  patterns  associated 
with  jet  streams,  apparent  intensity  of  tropical  storms,  and  so  forth. 
Great  impetus  was  given  to  this  method  with  the  launching  of  the  geo- 
stationary Applications  Technology  Satellites ,  namely  ATS- 1,  over  the 
Pacific  Ocean  equator  in  November  1967.    Time-lapse  photographs  are 
taken  with  the  Spin-Scan  Cloud  Camera  (SSCC)  approximately  every  25 
minutes  during  daylight  hours.    From  the  geostationary  altitude  of  approxi- 
mately 19,000    nmi,  the  camera  has  a  resolution  of  two  miles  at  zero 
nadir  angle,  and  views  a  geographical  region  from  about  50  S  to  50  N 
latitudinally  and  over  100  degrees  of  longitude  along  the  equator.    Conse- 
quently, excellent  daylight  video  coverage  of  cloud  systems  is  obtained 
over  the  oceanic  tropical  and  subtropical  areas  where  conventional  data 
are  normally  sparse.    Appendix  A  contains  ATS-i  photographs  taken  at 
approximately  1800  GMT  for  9-12  March  1971. 

Utilizing  methods  developed  primarily  at  the  Stanford  Research 
Institute  (SRI)   [19,  20,  21]  ,    NESS  personnel  produce  wind  estimates 
from  film  loops  of  ATS  time-lapse  photographs.    An  individual  scan  is 
superimposed  on  a  map  of  appropriate  scale  and  projection,  particular 
cloud  elements  are  selected  and  judged  to  be  either  upper  level  (near 
300  mb)  or  low  level  (near  850  mb) ,  and  marked  for  study.    These  cloud 
elements  are  then  followed  on  successive  frames  for  at  least  one  hour 
in  order  to  obtain  a  wind  estimate. 

An  example  of  the  operational  output  of  ATS  wind  data  from  NESS 
is  shown  in  Fig.   1.      These  data  are  incorporated  into  the   National 
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FIG.    1  300  MB.    MS  V'lNDS 


PROJECTION:   POLAR  STEREOGRAPHIC-TBUE  AT  60     NORTH    LATITUDE 
SCALE     160.000.000  ' 


FLEET  NUMERICAL  WEATHER  CENTRAL 
MONTEREY.  CALIFORNIA 


CHAAT  No  6B  Z 


Fig.   1.        0000  GMT  12  March  1971  ATS  winds 
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Meteorological  Center's    (NMC)    0000  GMT   Northern  Hemisphere  and 
tropical  analyses. 

C.        STATISTICALLY  DERIVED  500  mb  HEIGHTS  FROM  SATELLITE- 
OBSERVED  CLOUD  PATTERNS 

In  1968,  Nagle  and  Clark  [5,   18  ]  reported  on  a  quasi-objective 
method  of  incorporating  satellite-derived  cloud  pattern  information  into  a 
500  mb  numerical  analysis  scheme.    Subsequently,  Nagle  and  Hayden 
[10]    further  refined  this  technique  at  NESS,  developing  operational  pro- 
cedures used  by  the  NMC  quality  control  group  since  1969. 

The  NESS  technique  consists  of  the  following  procedures .    The  12- 
hour  500  mb  prognosis  from  1200  GMT,  which  forms  the  "first  guess" 
field  for  the  subsequent  0000  GMT  analysis,  is  subjected  to  ~  pattern- 
separation  program  developed    by  Holl  [6  3.      Two  of  the  several  fields 
produced  by  this  program  are  the  SR  or  residual  pattern,  which  may  be 
considered  a  space  mean  field  depicting  the  long-wave  pattern,  and  the 
SD  or  disturbance  field,  depicting  the  short-wave  features.    These  fields 
are  additive,  that  is 

Z    =    ZSR     +     2SD  (1) 

where  Equations  (1)  through  (5)  refer  to  500  mb. 

Assuming  the  SR  field  is  quasi-conservative  with  time  and  a  perfect 
prognosis,  the  SR  prognostic  value  accurately  specifies  the  current  (tg) 
SR  value.    Obtaining  the  current  500  mb  value  then  reduces  to  specifying 
the  current  SD  value  at  the  particular  location: 
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2  =      Z  +        Z  (?) 

tQ  SR  12  hr  prog  SDt  v  ' 

(t_12to  t0) 

Over  data-sparse  oceanic  regions,  conventional  observations  normally 
will  not  adequately  describe  the  fast-changing  SD  (tQ )  short-wavelength 
disturbances.    This  SD  value  is  obtained  by  statistically  linking  it  to 
characteristic  features  of  real-time,  satellite-observed  cloud  patterns: 

Z+      =    Zn_   ,_  .               +  Z0__.          (derived  from  current       (3) 

tn            SR  12  hr  prog  SDt          .       ,                 . 

/+         «.„  +    \  l0     cloud  pattern) 
(t_12  tot0) 

The  SD  field  approximates  that  of  the  geostrophic  relative  vorticity 
with  algebraic  sign  reversed.    Developing  cloud  systems  on  the  synoptic 
scale  are  related  to  the  vertical  motion  associated  with  the  positive  ad- 
vection  of  vorticity  in  the  upper  troposphere     [  9  ]   .      Thus ,  a  strong 
positive  correlation  has  been  obtained  between  negative  centers  in  the 

SD  field  and  centers  of  spiral  (vortical)  cloud  patterns    [5,   18,   10,  and 

2 
Mr.  F.  Bittner,     personal  communication,   1971  J    .    An  adjunct  of  this 

relationship  is  the  correlation  of  the  zero  line  of  the  SD  field  (i.e. ,  zero 

relative  vorticity)  with  cloud  patterns.    Thus,  superimposing  the  zero  SD 

line  upon  the  SR  field  immediately  yields  the  current  500  mb  height  from 

the  relationship 

Zt0  ZSR  12  hr  prog    +     ZSD_        ..  (4) 

.  .      -T  »  Zero  line . 

(t_12  to  tQ) 

=    z  (s) 

SR  12  hr  prog  K  ' 

(t_12totQ) 


2Mr.  Fred  Bittner  of  Project  FAMOS.  (This  is  the  popular  designator  for 
the  Naval  Weapons  Engineering  Support  Activity,  Iverson  Mall,  3737  Branch 
Ave.,  Room  307,  Hillcrest  Heights,  Maryland). 
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Fig.  2  gives  a  schematic  representation  of  many  of  these  correlated  rela- 
tionships between  the  satellite-observed  patterns  and  the  500  mb  SD  dis- 
turbance field. 

The  magnitude  of  the  SD  minimum  value  at  the  spiral  cloud  center 
(SGC)  is  obtained  from  the  solution  of  a  regression  equation  whose  pre- 
dictors are:    a)  latitude  of  the  SCC,    b)  maximum  diameter  of  the  SCC 
cloud  system,    c)  amplitude  of  the  SCC  cloud  system,  and  d)  the  Laplacian 
of  the  SR  field  at  the  SCC.    The  method  for  measuring  diameter  and  ampli- 
tude of  the  spiral  cloud  pattern  is  illustrated  in  Fig.  3.    An  alternate  re- 
gression equation  involving  only  the  latitude  of  the  SCC  and  the  Laplacian 
of  the  SR  field  at  the  SCC  also  has  been  developed  by  Nagle  and  Hayden 
for  situations  where  the  amplitude  and  diameter  of  the  cloud  pattern  are 
not  readily  obtained  from  the  satellite  picture. 

Both  of  these  regression  equations  are  shown  in  worksheet  form  in 
Fig .  4 .      The  predictor  coefficients  and  constants  in  the  worksheet  were 
given  to  the  authors  by  Mr.  Nagle  in  June  1971.     (These  coefficients  are 
derived  from  continuing  statistical  analyses  of  dependent  case  studies  of 
the  pattern- separated    components  of  the  NMC  numerical  500  mb  analysis 
and  spiral  cloud  systems  in  the  North  Atlantic  Ocean.) 

D.        SATELLITE  INFRARED  SPECTROMETER  (SIRS)  PRESSURE-TEMPERATURE 

SOUNDINGS 

Satellite  radiance  measurements  have  been  available  to  meteorolo- 
gists since  the  launching  of  the  first  TIROS  satellite  in  1960.    Although  the 
potential  of  these  observations  has  been  widely  studied,  their  incorporation 
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Fig.  2.      Schematic  representation  of  the  relationships  between  cloud 

patterns  and  features  of  the  500  rnb  relative  vorticity  field  [10]   . 


into  the  NMC  numerical  weather  analysis  scheme  Legal:  experimentally  in 
1969  with  the  conversion  of  Satellite  Infrared  Spectrometer  (SIRS)  radiation 
measurements  into  pressure-temperature  soundings     [16  ]  .    This  experi- 
mental work  has  since  become  a  real-time  effort  for  NESS  personnel  who 
now  routinely  process  these  SIRS  observations  for  use  by  the  quality  con- 
trol group  of  NMC  in  the  Northern  Hemisphere  analyses  for  both  0000  and 
1200  GMT  (Mr.  F.  Bittner,  personal  communication,   1971).    A  brief  sum- 
mary of  the  conversion  process  follows: 

The  SIRS  eight-channel  radiometer  observes  the  earth  and  its  atmo- 
sphere every  eight  seconds  along  the  north-south  track  of  the  NIMBUS 
satellite  from  an  orbital  altitude  of  approximately  600  nmi.    Its  field  of 
view  is  11.5  degrees,  or  a  "square"  of  the  earth's  surface  120  nmi  on  a 
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SINAP  WORKSHEET  (Fig.   4)  parameters: 
(§)  MAXDI  =    DIAMETER  (MM') 

©         CDAMP=    AMPLITUDE  (TT') 


Fig.  3.      Geometrical  construction  for  determining  the  diameter  and 
amplitude  of  a  cloud  system    [10  ]  . 
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Sq  as  the  grid  point  nearest  the 
SSC  specifies  the  LAPLACIAN 
heights: 

©    LAPSR=  (Si+  S2  +  S3  +  S4)  -   4  S0  = 
(E)    SSC     =    m 


r 

!2:- 


m 


m 


.S,       3"- 


m 


•  S. 


+/-  m 


°4- 

Sum  =  _ 

4xS0=. 


m 


m 


m 


REGRESSION  COMPUTATION  (SCC  METHOD) 

+ 
CONSTANT:  +21.5 


03 


x 
x 


(-  1.7): 

(-10.4): 
'  _  i     n\ , 

(-2.1): 
TOTALS:  +. 

SUM: 


ALTERNATE  COMPUTATION 
(Poor  Cloud  Dev.) 


+ 


CONSTANT 
@x    (-1.03); 


-53.5 


© 

©  SR  Value  at  SCC 

©  SD  Intensity  at  SCC 

(G)  Z500=SR+SD=  E    + 


+/-m 


F       = 


m 
m 
m 


© 

© 
© 


x    (-1.65): 
TOTALS:  +_ 

SUM: 


+/-m 


m 
m 
m 


Hydrostatic  HT.  Extrap.    :@Pq: 

©  Base  HT  (Table  IV)  *  :    
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*Tables  from  500  mb  height  extrapolation  method  from  Duthie ,  W.D.  , 
unpublished  lab  notes,  Meteorology  Dept.  ,  Naval  Postgraduate  School, 
Monterey,    Calif. 


Fig.  4a.   Worksheet  for  Nagle-Hayden  regression  equation  solution 
determining  500  mb  height  at  a  spiral  cloud  center  (SCC) 
using  satellite  pictures  and  the  SR500  field  on  a  1:30,000,000 
scale  (FNWC  chart  No.  4-P-l)  polar  stereographic  projection, 
true  at  6 ON. 
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Fig„  4b.     A  1/2  in  X  1/2  in  grid  over  an  FNWC  chart  No.  4-P-l  and 
compatible  with  worksheet  (Fig.  4a). 
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side.    These  readings  are  transmitted  to  CDA  stations  for  relay  to  the 
NESS  processing  center  at  Suitland,  Maryland.    At  NESS,   the  automated 
computer  processing  program  compares  three  successive  observations. 
From  these,    the  one  with  the  highest  sum  of  radiances  for  all  eight 
channels  is  selected.    Since  the  earth's  surface  is  usually  warmer  than 
cloud  tops,  this  procedure  tends  to  select  cloud-free  areas.     On  the 
average,  an  observation  is  selected  every  48  seconds,  or  every  three 
degrees  of  latitude,  along  the  track  of  the  satellite  subpoint.    Using  all 
available  radiance  readings  from  the  eight  channels  (channel  readings 
containing  excessive  noise  "spikes"  are  discarded),  temperature  profiles 
are  then  calculated  by  a  least-squares  multiple  regression  solution.    Best 
results  are,  of  course,  obtained  for  the  cloud-free  cases  when  all  eight 

channels  are  utilized.    The  solution  is  less  satisfactory  when  clouds  are 

3 
present  to  "contaminate"  channel  readings.       The  poorest  results  are 

obtained  for  high    cirrus    overcast      Computations  are  considered  unreli- 
able for  a  cirrus  coverage  greater  than  40  percent. 

Having  obtained  the  temperature  profile,  the  heights  of  constant- 
pressure  surfaces  are  computed  through  the  hydrostatic  relationship  by 

4 
using  the  850  mb  surface  as  the  known  reference  level.    The  final  result 

is  a  pressure-temperature  sounding  from  1000  mb  to  10  mb  in  a  format 

readily  usable  in  the  NMC  conventional  data  processing  program. 


^With  future  high  resolution  scanning  radiometers  ,  it  is  expected  that  tem- 
perature profiles  for  cloudy  conditions  will  be  as  accurate  as  the  present 
cloud-free  profiles.     Methods  are  also  being  derived  to  account  for  the 
influences  of  high  and  hot  terrain. 

^Depending  on  the  time  of  the  SIRS  observation,  interpolation  is  made 
between  the  preceding  850  mb  analysis  anc  the  850  mb  12-hour  prognosis 
to  obtain  the  most  accurate  850  mb  heignt. 
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At  NMC,  the  SIRS  soundings  are  subjected  to  the  same  accuracy 
checks  as  conventional  data;  that  is ,  any  soundings  differing  by  more  than 
200  m  when  compared  with  the  500  mb  and  300  mb  "first  guess"  heights 
are  rejected  in  their  entirety.    Such  rejected  soundings  are  subjectively 
reviewed  by  a  quality  control  group  which  may  re-insert  the  sounding  if  it 
is  deemed  accurate.    All  accepted  soundings  are  put  into  the  objective 
analysis  program.    The  SIRS  soundings  are  given  the  same  "weight"  or 
significance  as  the  conventional  soundings. 

NMC  performs  both  an  operational  upper  air  analysis  (at  approxi- 
mately 0000/1200  GMT  +3.5  hours)  and  an  update  analysis  (0000/1200 
GMT  +  10  hours).    The  latter  utilizes  all  the  synoptic  or  near  synoptic 
data  received  since  0000/1200  GMT.    Thus,  the  SIRS  soundings, which  are 
being  continuously  processed,  are  usually  entered  into  both  analyses. 
Presently,  the  processing  time  delay  at  NESS  dictates  that  only  those 
SIRS  observations  taken  in  the  interval  from  one  to  six  hours  prior  to  synop- 
tic time  (0000/1200  GMT)  can  be  processed  in  time  for  incorporation  into 
the  operational  analysis.    The  number  of  SIRS  soundings  entering  the 
operational  analysis  is  usually  less  than  30;  whereas  for  the  update 
version,  the  total  (including  those  in  the  operational  analyses)  averages 
about  60  soundings.    The  comparable  average  figures  for  conventional 
radiosonde  data  are  500  and  575,  respectively  (March  1971  figures). 

An  example  of  the  NMC  update  500  mb  SIRS  coverage  for  0000  GMT 
12  March  1971  is  shown  in  Fig.  5.    The  lack  of  data  at  lower  latitudes  is 
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Fig.  5.      0000  GMT  12  March  1971   500  mb  SIRS  data  with  the  analyzed 
FNWC  heights  in  parentheses. 
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partly  a  function  of  the  wider  spacing  of  orbital  paths  near  the  equator. 
As  the  orbital  paths  converge  in  the  polar  areas,  the  coverage  increases 
accordingly. 


Future  satellite  spectrometers  will  also  scan  in  a  direction  lateral  to 
their  orbital  paths,  providing  data  to  fill  in  these  wide  gaps  in  the  lower 
latitudes. 
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III.      DATA  SOURCES  AND  SELECTION  OF  STUDY  PERIOD 

The  proposed  interaction  with  the  FNWC  analysis  and  prediction 
scheme  hinged  on  their  retention  (on  magnetic  tape)  of  the  conventional 
synoptic  data  so  that  a  re-analysis  could  be  performed  at  a  subsequent 
date  using  bogus  reports  fabricated  from  the  satellite-derived  information. 
FNWC  agreed  to  save  their  update  synoptic  data  during  the  month  of 
March  1971,  thereby  establishing  the  month  of  the  study  period. 

The  National  Environmental  Satellite  Service  (NESS)  was  contacted 
by  Mr.  F.  Bittner  of  Project  FAMOS  for  use  of  their  ATS  film  loops  and 
photographs,  AVCS  video,  and  ITOS  nighttime  IR  Northern  Hemisphere 
mosaics .     However,   Lo  ensure  a  ye^e  satellite  date:  bc>se  froir.  v/hich  i-o 
work,  in  the  eventuality  of  the  AVCS  or  ATS  data  from  NESS  not  being 
available,  a  large  number  of  satellite  data  sources  were  requested  to  pro- 
vide their  satellite  material. 

All  shore-based  Naval  Weather  Service  Command  activities  in  the 
North  Pacific  Ocean  having  satellite  readout  equipment  were  contacted 
for  possible  use  of  their  gridded  pictures.     Lt.  Cmdr.  D.  Edgren  of 
NWSED,  Moffett  Field,  Calif.  ,  contacted  the  Air  Force  Satellite  Center, 
Sunnyvale,  Calif.,  for  their  APT  readouts.    Cmdr.  W.  Willingham  of 
FWF  Suitland,  Md .  ,  contacted  the  NMC  for  the  operational  SIRS  sound- 
ings and  nephanalyses  .    The  response  of  all  activities  was   gratifying 
and  ensured  that  an  excellent  collection  of  satellite  pictures  was  available. 
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The  five  upper  air  synoptic  times  from  0000  GMT  10  March  through 
0000  GMT  12  March  1971  were  selected  for  the  study  period  for  two 
reasons: 

1)  neither  satellite-derived  bogus  data  from  FWF  Suitland 
nor  SIRS  data  from  NMC  (via  FWF  Suitland)  entered  the 
FNWC  analyses  during  this  period. 

2)  substantial  short  period  weather  pattern  changes  were 
observed  in  the  North  Pacific  Ocean. 

The  first  reason  ensured  that  the  FNWC  objective  analyses  were 
essentially  void  of  subjective  or  satellite-influenced  changes.    To  further 
this  concept,  the  real-time  surface  cyclone  boguses  incorporated  in  the 
operational  FNWC  analyses  from  the  FWC  Pearl  Harbor  were  deleted  from 
the  magnetic  tape  data  retention.    Thus,  in  the  later  re-running  of  this 
data,  both  an  UNMOD  ANAL  containing  only  conventional  data,  and  a 
MOD  ANAL  containing  conventional    and  the  authors'    bogus  insertions 
were  generated . 

The  second  reason  was  based  on  the  premise  that  the  weather  pattern 
changes  noted  on  the  synoptic  maps  would  be  noticeable  in  the  satellite 
data. 

Table  I  lists  the  satellite  and  conventional  data  sources  used  in  this 
research.    The  authors  are  indebted  to  each  activity  and  the  individuals 
in  these  organizations  for  their  cooperation  in  making  this  data  available 
to  the  Department  of  Meteorology,  Naval  Postgraduate  School,  Monterey, 
California. 
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TABLE   I. 
Sources  of  conventional  and  satellite  data 

A.  Satellite  Data 

1.         Video  Pictures  (ITOS-1,  ESSA  VIII)  : 
FWC  Alameda 
FWC  Pearl  Harbor 
FWC  Guam 
FWF  Kodiak 
FWF  Yokusuka 
FWF  Sangley  Point 
NPS  Monterey 
PMR  Pt.  Mugu 
NWSED  Moffett  Field 
tvtw^eD  Midwav  T^^^^ 
NWSED  Adak  Isc 
AFSC  Sunnyvale 
NESS,  Suitland  (Mosaics) 

*2.         IR  Nighttime  Data  (ITOS-1)  : 
NESS,  Suitland  (Mosaics) 

*30         ATS-1  Wind  Computations  and  Film  Loops  : 
NESS,  Suitland 

*4.         SIRS  Soundings: 

NESS,  Suitland 

NWSED  Asheville  (repository  for  NESS  SIRS  soundings) 

*5.  Nephanalyses: 

NESS,  Suitland  (also  available  from  FOFAX  and 
NAFAX  circuits) 

*6.         Pacific  Operational  Bogus  Data  and  SINAP  Maps  : 

NMC,  Suitland 
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B.  Conventional  Data 

1.         Surface  and  Upper  Air  Synoptic  Data: 

FNWC  Monterey 
NMC,  Suitland 

2  .         Surface  Analyses  ; 

FNWC  Monterey 

FWC  Alameda 

FWC  Pearl  Harbor 

NMC,  Suitland  (FOFAX  and  NAFAX  circuits) 

3  .         Upper  Air  Analyses  : 

FNWC  Monterey 
FWC  Pearl  Harbor  (tropical) 
NMC,  Suitland  (FOFAX  and  NAFAX  circuits) 
4.         Pattern  Separation  Analyses  (SR  and  SD  Contour  Fields) 
FNWC  Monterey 

*Collected  and  forwarded  by  Project  FAMOS 
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IV.      PROCEDURES  AND  RESTRICTIONS 

Quasi-operational  ground  rules  were  established  for  incorporating 
the  satellite-derived  data  into  the  data  base  of  the  FNWC  numerical 
analysis  program.    The  primary  rule  was  that  all  procedures  be  designed 
to  improve  the  operational  FNWC  analyses  (available  at  0000/1200  GMT 
+   4.5  hours)  in  time  for  these  improvements  to  enter  the  update  analysis 
run  (approximately  0000/1200  GMT  +  10  hours).    This  implied  that  the 
following  intermediate  (0600/1800  GMT)  surface  analysis  (produced  at 
0600/1800  GMT  +   2  hours)  was  available  to  indicate  developing  trends; 
and  the  operational  0000/1200  GMT  500  mb  SR  pattern  could  be  used  in- 
stead of  the  previous  12-hour  SR  prognosis  as  used  in  the  NESS  version 
of  the  Nagle-Hayden  technique. 

Another  important  rule  was  that  only  satellite  pictures  in  the  synop- 
tic time  frame  could  be  used.     For  example,  APT  photographs  taken  after 
the  1000  GMT  cut-off  time  for  the  0000  GMT  analysis  were  not  used  in 
modifying  the  original  analysis.    In  the  actual  research  this  rule  was  not 
a  hindrance  as  most  Pacific  video  satellite  pictures  are  taken  between 
2000  GMT  and  0200  GMT;  thus,  they  are  naturally  restricted  to  usage  in 
the  0000  GMT  analysis.    The  only  views  available  for  the  1200  GMT 
synoptic  period  were  the  ITOS  IR  nighttime  mosaics  taken  between  1000 
GMT  and  1600  GMT  in  the  Pacific  Ocean  region. 
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The  third  operational  rule  was  that  conventional  observations  were 
considered  reliable  and  inviolate  unless  clearly  demonstrated  otherwise. 

Certain  restrictions  were  placed  on  the  re-analysis  effort  owing  to 
the  current  (July  1971)  programming  practices  at  FNWC;  namely,  bogus 
inputs  were  permitted  only  at  sea   level  (pressure  and  wind)  and  at  500  mb 
(height,  temperature  and  wind).     Furthermore,  the  total  number  of  bogus 
inputs  was  limited  to  100  reports  at  sea  level  and  to  50  reports  at  500  mb. 
This  limit  of  fifty  boguses  at  500  mb  proved  a  handicap  in  achieving  the 
desired  re-analysis  at  this  level. 

Within  these  rules  and  restrictions,  the  following  procedures  were 
used  to  modify  the  original  numerical  analyses  utilizing  the  weather- 
satellite  data  up  to  the  equivalent  real-time  data  cut-off  point  for  each 
re-analysis: 

1)  The  operational  0000/1200  GMT  surface,  500  mb,  500  mb  SR, 
300  mb,  and  thickness  (500  mb  -  1000  mb)  fields  were  considered  the 
first  guess  fields . 

2)  SIRS  500  mb  height  data  and  ATS-1   300  mb  wind  data  were 
plotted  on  the  corresponding  analyses,  and  significant  differences  noted. 

3)  All  available  satellite  cloud  pattern  pictures  were  reviewed 
to  gain  a  familiarity  between  cloud  patterns  and  synoptic  systems. 


The  FNWC  program  has  since  been  rewritten  to  accommodate  100 
boguses  at  sea  level,  500  mb    and   300   mb   (Mr.    T.    Hesse, 
FNWC,  personal  communication,  August  1971). 
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4)  Using  the  500  mb  SR  analysis,  spiral  cloud  centers  and  zero 
vorticity  (SD  zero-line)  lines  noted  on  the  satellite  views  were  plotted 
and  the  Nagle-Hayden  statistical  computations  performed. 

5)  The  FNWC  surface  and  500  mb  analyses  were  subjectively 
changed  to  conform  to  the  major  features  derived  from  the  satellite  data. 

6)  The  next  0600/1800  GMT  surface  analysis  was  used  to  im- 
prove the  sea-level  modification. 

7)  A  1000  mb  field  was  constructed  from  the  modified  surface 
charts,  then  a  new  thickness  was  manually  generated. 

8)  The  subjective  thickness  was  compared  to  the  numerical  12- 
hour  history  thickness.    Any  changes  in  the  normally  conservative  thick- 
ness patterns  were  reconciled,  with  necessary  adjustments  made  to  the 
appropriate  thickness  pattern.     From  the  changed  thickness  pattern,  the 
surface  and  500  mb  analyses  were  adjusted  to  reflect  the  change  and  to 
ensure  continuity  of  the  pressure  systems. 

9)  Bogus  reports  (see  Appendix  C)  were  prepared  for  the  purpose 
of  objectively  duplicating  the  manually  modified  surface  and  500  mb 
charts . 

10)  Analyses  were  then  made  of  the  surface  and  mandatory  upper 
air  levels  using  only  the  conventional  data  first  (UNMOD  ANALS),  and 
then  using  both  the  conventional  and  bogus  data  (MOD  ANALS). 

11)  The  re-analyses  were  examined.    If  cloud-related  synoptic 

features  were  not  sufficiently  delineated,  new  bogus  data  were  fabricated,  and 

a  re-cycling  performed  until  the  modified  fields  appeared  satisfactory 
(experience  would  tend  to  eliminate  the  re-cycling). 
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12)  The  primitive  equation  (PE)  prognosis  was  run  to  48  hours  on 
both  the  unmodified  and  satellite-modified  analysis  fields,  generating 
UNMOD  and  MOD  PROGS. 

13)  Finally,  error  fields  for  the  surface  and  500  mb  levels  were 
generated  for  both  sets  of  PE  prognostics. 
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V.       CASE  STUDIES  OF  MODIFICATIONS  TO  ANALYSES 

Synoptic  times  0000  GMT  10  March  through  1200  GMT  12  March 
1971  were  studied  in  detail,  with  modifications  made  to  the  respective   sea 
level  pressure  and  500  mb  fields.    The  dates  of  1200  GMT  10  March  and 
0000  GMT  11  March  1971  were  selected  as  initial  times  from  which  the 
current  (July  1971)  FNWC  primitive  equation  (PE)  five-layer  prognostic 
model  was  run  to  48  hours.    Several  changes  had  been  made  to  the  prog- 
nostic program  since  March  1971;  therefore,  to  make  a  valid  comparison 
of  the  UNMOD  PROG  and  the  MOD  PROG,  the  UNMOD  PROG  was  re-run. 

Figures  showing  analyses  follow  immediately  after  the  detailed 
description  of  the  ana-Lysis  modifications  made  in  each  c:acc  study.    Tht; 
vortices  are  labeled  alphabetically  from  Eastern  Pacific  to  Western  Pacific 
as  they  made  their  appearance  in  time.    In  each  case,  the  referenced 
positions  are  the  500  mb  position  unless  otherwise  stated.    The  dashed 
contours  of  the  analyses  are  isolines  of  the  UNMOD  ANAL;  they  are  super- 
imposed on  the  neoteric  MOD  ANAL  isolines  to  clearly  illustrate  the  modi- 
fications effected  by  the  inclusion  of  explicit  satellite  data  and  the  implicit 
"values"  of  the  satellite  observations. 

A.         0000  GMT  10  MARCH  1971  MODIFICATIONS  TO  ANALYSES 

The  0000  GMT  surface  maps  of  the  Pacific  Ocean  are  based  on  a 
maximum  number  of  ship  reports  because  it  is  near  midday  over  this  region 
and  therefore  fits  into  the  routine  broadcast  day.    Video  satellite  pictures 
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are  available  at  this  time  and  are  the  basis  of  the  following  modifications 
to  the  sea  level  and  500  mb  MOD  ANALS  (Figs.  6,  7  and  8): 

1.  Vortex  A 

Vortex  A  at  62N  134W  is  over  British  Columbia  and  is  in  the 
dissipation  stage  based  on  the  fragmented  appearance  of  its  spiral  cloud 
pattern . 

2.  Vortex  B 

Vortex  B,  near  5 IN    156W,  is  a  weak  trough  on  the  UNMOD 
ANAL.     In  the  video  pictures,  this  vortex  appears  to  be  a  classical  mature 
cyclonic  vortex  associated  with  a  major  polar  frontal  band.    Consequently, 
at  500  mb  a  low  with  a  central  value  of  5000  m  was  introduced.    Its  asso- 
ciated trough  was  repositioned  westward  on  the  basis  of  the  sharp  dissi- 
pation of  middle  and  high  clouds  in  the  frontal  band  and  the  curvature  of 
the  cloud  lines . 

3.  Vortex  C 

Vortex  C  at  52N    179E,    with  a  central  value  of  5072  m  on  the 
UNMOD  ANAL,  was  relocated  to  50N    176E   based  on  the  satellite  picture 
of  the  center  and  lowered  to  5000  m  based  on  the  computational  values  ob- 
tained.   This  is  a  mature  cold  vortex.    Its  trough  was  deepened  to  the 
south  and  southwest  and  the  cyclonic  curvature  accentuated  on  the  basis 
of  the  cloud  band  curvatures  . 

4.  Vortex  D 

Vortex  D,  over  the  north  tip  of  Honshu,  is  well  defined  by 
Japanese  radiosonde  reports.      East  of  this  vortex,  the  ridge  along  155E 
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was  increased  slightly  based  on  the    anticyclonic  cloud   curvature  and 
the  clear  area  on  the  east  side  of  the  ridge. 

5.  Vortex  F 

A  comma  cloud  (PVA  Max)  at  4  IN    177W  located  northwest  of 
a  major  bulge  in  the  frontal  band,  presages    future  development  of  a  wave 
cyclone  (vortex  F).    The  trough  and  downstream  ridge  were  enhanced 
slightly  on  the  basis  of  the  SD  zero-line  delineation  and  thickness  pattern 
considerations. 

6.  Vortex  G 

Vortex  G,  over  the  northeast  tip  of  Siberia,  is  defined  by 
conventional  radiosonde  data;  the  video  pictures  contained  too  much  snow 
and  ice  reflection  to  permit  observation  of  the  cloud  systems. 

B.  1200  GMT  10  MARCH  1971  MODIFICATIONS  TO  ANALYSES 

The  1200  GMT  sea  level  pressure  map  of  the  Pacific  Ocean  normally 
has  few  ship  observations  because  it  is  near  midnight  (local  time)  over 
this  region,  and  therefore  most  ships  either  do  not  record  or  do  not  trans- 
mit observations  at  12  00  GMT  in  the  Pacific  Ocean.    Based  primarily  on 
the  photo  mosaic  of  ITOS    1   IR  nighttime  viewing,  the  following  modifica- 
tions were  made  to  the  sea  level  and  500  mb  MOD  ANALS  (Figs.  9,  10  and 

11): 

1.         Vortex  B 

Vortex  B  has  passed  through  the  mature  stage  based  on  the 
fracturing  of  the  dry  tongue  and  the  gray  tone  of  the  encircling  clouds  in 
the  IR  data.      However,  it  is   still  a  well-defined   closed   system.    The 

39 


original  FNWC  analysis  rejected  Ship  Papa's   500  mb  wind  of  335  deg/38 
kt,  as  it  failed  the  consistency  check  between  the  700  mb  wind  of  245  deg/ 
46  kt  and  the  400  mb  wind  of  240  deg/36  kt.    In  the  modification,  the 
500  mb  wind  was  changed  to  235  deg  in  consonance  with  the  cloud  pattern 
and  the  speed  unchanged.    Vortex  B  was  positioned  northwest  of  Ship  Papa 
with  an  estimated  central  height  of  4990  m.    However,  due  to  the  difficulty 
encountered  in  bogusing  (see  discussion  in  Appendix  C),  the  MOD  ANAL 
value  in  the  numerical  analysis  version   is     5047  m.    To  maintain  vertical 
consistency  and  continuity,  the  sea-level  low  was  positioned  beneath  the 
500  mb  low  and  deepened  from  980  mb  to  971  mb.    The  500  mb  trough  ex- 
tending southeast  from  vortex  B  was  deepened  slightly  in  the  vicinity  of 
the  jet  maximum  along  the  north  boundary  of  the  frontal  band  near  42N  14 0W, 

2.  Vortex  C 

Vortex  C,  south  of  Adak,  Alaska,  is  probably  nearing  the  dis- 
sipating stage  as  the  brightest  clouds  (coldest  temperatures)  completely 
surround  the  center,  and  the  cloud  lines  are  generally  fragmented  and  not 
bright,  indicating  stable  mid-  or  low-level  stratocumulus .    Winds  are 
missing  from  Adak's  report  (station  70454),  but  the  clouds  detected  by  the 
IR  sensor  indicate  weak  easterly  flow.    The  500  mb  analysis  was  modified 
to  show  this  closed  circulation. 

3.  Vortex  D 

Vortex  D,  east  of  Hokkaido,    has  not  intensified  based  on  both 
conventional  data  and  cloud  patterns .    The  bright  clouds  nearly  surrounding 
the  center  indicate  this  small  system  is  in  the  mature  stage  „ 
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4.  Vortex  E 

A  large,  bright  cloud  mass  southeast  of  vortex  D,  near  40N 
158E,  has  separated  from  this  vortex  and  is  in  the  initial  stage  of  surface 
cyclone  development.    Vortex  E  is  not  as  well  organized  as  vortex  F  at 
this  time,  but  the  photo  of  IR  data  clearly  shows  that  this  is  the  next 
unstable  system  upstream  from  vortex  F.    Consequently,  at  500  mb  the 
trough  on  the  west  side  of  the  cloud  mass  was  deepened  and  the  ridge  on 
the  east  side  (near  162E)  was  enhanced  by  using  the  SD  zero- line  to  de- 
rive 500  mb  heights.    However,  the  numerical  MOD  ANAL  did  not  portray 
the  short-wave  character  intended  for  this  system. 

5.  Vortex  F 

The  comma  cloud  noted  in  the  0000  GMT  10  March  discussion 
has  apparently  joined  with  the  convex  bulge  in  the  polar  front  near  42N 
16 0W.    This  is  between  the  two  major  vortices,  B  and  C,  in  an  area  of 
maximum  baroclinicity.    The  photo  of  IR  values  shows  this  cloud  mass  to 
be  the  whitest  region  of  cloudiness  in  this  portion  of  the  ocean.    A  frontal 
wave  which  exhibits  an  area  of  cold  temperatures  (white)  in  the  IR  mode, 
suggesting  strong  ascending  motion  and  falling  sea  level  pressures,  is 
indicative  of  a  developing  surface  cyclone  and  is  typical  of  the  second 
stage  of  a  rapidly  developing  system    [15]   .      Cyclonic  circulation  is 
weakly  evident,  though  the  dry  tongue  is  visible.    This  cloud  pattern  is 
nearly  identical  to  the  CP20  example  of  Nagle  and  Clark,  which  they  de- 
scribe as  indicating  a  rapidly  deepening  system.    Thus  the  pattern  of  IR 
values  of  vortex  F  implies  a  strong  short-wave  trough  at  500  mb  with  a 

closed  surface  circulation. 
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It  is  interesting  to  note  that  the  sea  level  pressure  could  not  be 
properly  ascertained  in  this  area  from  conventional  ship  reports,  the 
nearest  reporting  ship  being  over  500  nmi  away.    Thus  the  FNWC  surface 
analysis  has   no  indication  of  this  developing  system.    Even  the  1800  GMT 
10  March  surface  analysis  continued  to  show  fairly  zonal  flow  over  the 
region.    In  fact,  a  key  ship  report  (call  sign  JIVA)  near  50N  150W  reported 
a  pressure  of  984.5  mb  and  south-southwest  wind  of  40  kt.    Its  pressure 
was  interpreted  to  belong  in  the  pattern  of  vortex  B  to  the  northeast, 
though  the  wind  indicated  cyclonic  development  to  the  west  of  the  ship. 
Unfortunately,  no  other  ship  observations  were  available  so  the  FNWC 
duty  officer  attached  a  note  to  this  report,    "...  in  without  wind,"    indi- 
cating that  he  considered  the  pressure  correct,  but  the  wind  erroneous. 
Not  until  0000  GMT  11  March  did  six  ships  near  the  storm  send  in  their 
reports  thereby  pinpointing  the  storm's  position  and  intensity  (975  mb). 
The  National  Weather  Service  maps  were  equally  vague  in  defining  this 
system  until  0000  GMT  11  March. 

The  Nagle-Hayden  technique  produced  500  mb  values  of  5320  m  for 
the  spiral  cloud  method,  and  5370  m  for  the  alternate  method.    Extrapo- 
lating back  in  time  from  a  re-analysis  of  the  1800  GMT  surface  analysis 
indicates  that  the  sea  level  pressure  was  below  996  mb  at  1200  GMT, 

thus  in  hydrostatic  agreement  with  the  final  subjective  500  mb  height  of 

7 
5340  m.         The  original  analysis  value  at  this  point  was  5460  m.    In  the 


7 

Note  that  the  1800  GMT  surface  analysis  is  available  prior  to  the 

update  of  the  1200  GMT  analysis  at  FNWC. 
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modification,  no  attempt  was  made  to  close  off  a  500  mb  cyclonic  circu- 
lation.   The  500  mb  trough  was  re-positioned  on  the  upstream  edge  of  the 
cloud  pattern  associated  with  the  developing  vortex  and  deepened  accord- 
ing to  the  height  computation  of  5340  m.    The  sea-level  low  was  positioned 
by  the  apparent  circulation  center  in  the  cloud  patterns  and  deepened  to  a 
closed  995  mb.     (Further  discussion  of  this  vortex  is  contained  in 
Section  VI.) 

6.         Vortex  G 

Vortex  G  at  62N    166E  again  is  not  distinguishable  on  the 
photo  of  IR  data,  probably  due  to  very  low  clouds  or  an  isothermal  atmo- 
sphere, but  was  adequately  covered  by  Russian  radiosonde  stations. 

C.         0000  GMT  II  MARCH  1971  MODIFICATIONS  TO  ANALYSIS 

An  excellent  suite  of  video  pictures  for  this  period  (1900  to  0200 
GMT)  was  obtained  from  Naval  Weather  Service  facilities  and  NESS.    The 
NWSED  Midway  Island  APT,  and  the  NESS  ATS-1  film  loops  and  enlarged 
photographs  were  particularly  useful  in  showing  the  rapid  development  of 
vortex  F. 

The  following  modifications  were  made  to  the  surface  and  500  mb 
analyses  (Figs.   12,   13  and  14): 
1.         Vortex  B 

Vortex  B  in  the  Gulf  of  Alaska  is  considered  to  be  in  a  late 
mature  stage.    Its  frontal  occlusion  has  separated  from  the  vortex  center, 
the  dry  tongue  nearly  encircles  the  center,  and  the  feeder  bands  are  short 
and  diffuse. 
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Ship  Papa's  (50N  145W)  wind  direction  at  500  mb  is  reported  as 
265  deg/54  kt,  but  the  cloud  streaks  on  the  video  pictures  indicate  west- 
northwesterly  flow  at  mid-troposphere  levels.    An  examination  of  the 
sounding  showed  a  constant  wind  direction  of  265  deg  from  850  mb  to 
100  mb,  which  is  not  likely  to  occur.    It  was  decided  to  insert  a  bogus 
west-northwesterly  wind  of  290  deg  in  the  vicinity  of  Ship  Papa.    Since 
the  conventional  report  could  not  be  modified  by  the  bogus  method  without 
rejecting  the  entire  sounding,  several  bogus  winds  at  500  mb  were  in- 
jected to  "override"   (see  Appendix  C)  the  Papa  report,  however,  the 
numerical  MOD  ANAL  achieved  only  a  westerly  flow  of  275  deg. 

The  desired  500  mb  central  value  for  Vortex  B  was  4980  m.    Deliber- 
ate "overbogus"         (see  Appendix  C)  of  four  surrounding  points  brought 
the  central  value  down  to  only  5061  m  but  achieved  the  desired  cyclonic 
circulation.    The  associated  trough  was  deepened  to  the  southeast  of  the 
vortex  as  indicated  by  the  cloud  curvatures  and  SD  zero-line  computations. 

2.         Vortex  C 

The  circulation  associated  with  vortex  C  is  well  defined  and 
covers  a  large  area  as  it  continues  to  drift  eastward  with  minor  intensifi- 
cation.   The  500  mb  center  was  closed  off  and  the  trough  deepened,  con- 
forming to  the  cyclonic  curvature  of  the  cloud  lines  and  to  the  change  in 
cloud  type  (open  cells  on  the  west  side  and  enhanced  cloudiness  on  the 
east  side  of  the  trough  line). 
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3.  Vortex  D 

Vortex  D  has  moved  slowly  eastward  and  has  weakened, 
judging  from  the  large  dry  tongue  which  completely  encircles  the  center. 

4.  Vortex  E 

Vortex  E  has  formed  between  vortex  D  and  the  large  bright 
frontal  cloud  mass  noted  in  the  1200  GMT  10  March  discussion.    As  a 
result  of  the  interaction  between  the  comma  cloud  (PVA  Max)  and  the  en- 
hanced convective  activity  along  the  frontal  zone,  wave  development  on 
the  front  is  indicated  in  the  vicinity  of  40N  160E.    The  surface  cyclone 
was  deepened  and  separated  from  the  surface  cyclone  associated  with 
vortex  D  to  the  northwest  to  indicate  its  interpreted  development.    At 
500  mb,  the  trough  was  deepened  nearly  70  m  and  the  downstream  ridge 
was  heightened  40  m.    Eight  500  mb  bogus  wind  reports  were  entered  on 
the  east  side  of  this  trough  to  attempt  depiction  of  the  strong  southwester- 
ly flow  in  this  region,  but  again  the  numerical  MOD  ANAL  could  not  in- 
corporate this  shortwave  intensification  as  was  desired. 

5.  Vortex  F 

As  noted  in  the  1200  GMT  10  March  discussion,  vortex  F,  in 
the  vicinity  of  48N  155W,  has  undergone  rapid  intensification  during  the 
night.    The  first  daylight  pass  of  ESSA  VIII  (orbit  #10233  at  1945  GMT 
10  March)  clearly  shows  the  intense  nature  of  this  small  system.    The 
cloud  pattern  corresponds  perfectly  with  the  Nagle  and  Clark  CP20  or 
"hammerhead"  designation    [5  ]  .       Vortex  F  must  be  considered  fully 
occluded  at  this  time,  dictating  a  500  mb  closed  low  slightly  northwest 
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of  the  sea-level  low,  flanked  by  strong  ridges  east  and  west,  south  of  a 
vigorous  trough  in  the  cold  air  along  155W.    Since  the  conventional 
analysis  had  a  slight  ridge  in  this  vicinity  and  only  minor  troughing  in 
the  nearby  systems  B  and   C ,  considerable  modification  was  required  to 
fit  the  dynamics  and  flow  pattern  of  vortex  F.    The  500  mb  heights  were 
deepened  as  much  as  300  m  in  the  trough.    The  500  mb  center  was  modi- 
fied to  show  a  closed  circulation  and  a  center  value  of  5090  m.      The 
ridges  upstream  and  downstream  were  heightened  by  using  the  SD  zero- 
line  computations . 

6.         Vortex  G 

Though  Vortex  G  is  extremely  small,  the  video  photo  shows 
it  to  be  a  mature,  intense  system.    Vortex  G  is  moving  southward,  enter- 
ing the  Bering  Sea  east  of  the  Kamchatka  Peninsula  in  conformity  with  the 
spiral  cloud  pattern.      The  Russian  data  network  measured  this  system 
very  accurately  and  no  modifications  were  required. 

D.         1200  GMT  11  MARCH  1971  MODIFICATIONS  TO  ANALYSES 

The  1200  GMT  11  March  1971  synoptic  maps  again  illustrate  the 
lack  of  conventional  data  in  the  North  Pacific.  Three  station  vessels, 
November,  Papa,  and  Victor,  were  the  only  rawinsonde  ships  reporting 
over  the  entire  oceanic  expanse.  Marcus,  Wake,  and  Midway  Islands 
were  the  only  island  stations  in  strategic  locations;  however,  the  500  mb 
winds  were  missing  in  the  reports  from  both  Midway  and  Wake.  Surface 
synoptic  ship  reports  were  equally  sparse.    In  the  eastern  North  Pacific 
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Ocean  between  180  longitude  and  the  west  coast    of  North  America,  eleven 
ship  reports  were  plotted  on  the  NWS  preliminary  surface  analysis,  and 
eight  ship  reports  were  plotted  on  the  FNWC  operational  surface  analysis. 

The  photo  mosaic  of  ITOS-1  IR  nighttime  viewing  was  the  primary 
source  for  the  following  changes  to  the  surface  and  500  mb  analyses 
(Figs.   15,    16  and  17): 

1.  Vortex  B 

Vortex  B  moved  northwestward  into  the  northern  Gulf  of  Alaska 
and  is  well  delineated  by  conventional  rawinsonde  observations  at  Kodiak 
(70350)  and  Anchorage  (70273).  It  is  poorly  defined  in  the  IR  data  values 
and  is  probably  being  absorbed  into  the  major  low  pressure  system  over 
northern  Alaska.  Its  trough  extends  southeastward  and  affects  the  circu- 
lation over  Ship  Papa  which  reported  a  500  mb  wind  of  235  deg/14  kt.  No 
modifications  were  required  for  vortex  B. 

2.  Vortex  C 

The  photo  mosaic  was  badly  blurred  in  the  vicinity  of  vortex 
C.    The  trough  position  was  discernible  from  the  cloud  patterns,  but  the 
center  was  not.    Relying  on  continuity,  vortex  C  was  considered  a  500  mb 
cyclone  with  a  central  height  of  5000  m. 

3.  Vortex  D 

Vortex  D,  east  of  Hokkaido,  was  completely  filled  with  gray 
clouds  (low  stratocumulus) ,  indicating  dissipation.    The  good  data  cover- 
age of  the  surface  specified  a  1016  mb  closed  sea-level  low.     The  surface 
analysis  was  not  changed.      Based  on  the  observed   cloud   lines, 
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indicating  primarily  the  orientation  of  thickness  contours,  the  500  mb 
analysis  was  modified  slightly  to  enhance  the  cyclonic  flow  south  and 
west  of  vortex  D., 

4.  Vortex  E 

Vortex  E  is  the  dominant  mid-Pacific  feature.    The  center  is 
faintly  discernible  on  the  mosaics  with  a  large,  extremely  wide  dry  tongue 
extending  northeast.    The  stage  of  development  is  early  maturity  (occluding) 
or  early  CP  20     [5]   .      In  the  area  of  cold  air  advection  (the  dry  tongue), 
large  amounts  of  gray  clouds  are  present,  suggesting  low  cloud  tops, 
stable  conditions,  and  a  small  air/sea  temperature  difference.    This 
would  indicate  that  deepening  is  proceeding  at  a  slow  rate    [13]  . 

Vortex  E  is  closed  off  with  a  5160  m  central  value  at  500  mb. 
The  surface  low  is  positioned  to  the  southeast  and  deepened  to  1004  mb. 
The  ridge  along  178E  has  begun  to  move  in  such  a  way  as  to  increase  the 
distance  between  it  and  vortex  E.    This  ridge  was  enhanced  to  accommo- 
date greater  anticyclonic  flow  apparent  from  the  cloud  patterns . 

5.  Vortex   F 

Vortex  F  is  the  dominant  East  Pacific  synoptic  feature  at  this 
time.    From  the  photo  of  IR  values  vortex  F  can  be  classified  as  a  CP30 
[5l    or  a  mature  storm  with  a  secondary  frontal  wave    [7  1.      Vortex  F  is 
approximately  350  nmi  southeast  of  Ship  Papa  with  the  sea  level  pressure 
lowered  to  964  mb  based  on  the  continuing  development.    The  500  mb 
center  was  lowered  to  5040  m  based  on  the  quasi-objective  computations. 
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6 .  Vortex  G 

Vortex  G  has  moved  further  into  the  Bering  Sea.    The  added 
moisture  supply  has  clearly  enhanced  the  cloud  development  seen  in  the 
ITOS  IR  depiction.     Minor  modifications  were  made  at  500  mb  to  increase 
the  cyclonic  flow  in  the  trough  extending  southwest  of  vortex  G. 

7.  Vortex  H 

An  incipient  secondary  frontal  wave  development  (Vortex  H) 
can  be  anticipated  because  of  the  enhanced  comma  cloud  (PVA  Max)  near 
40N  150W  on  the  IR  mosaic. 

E.  0000  GMT  12  MARCH  1971  MODIFICATIONS  TO  ANALYSES 

The  excellent  video  coverage  at  this  synoptic  time  permitted  a  de- 
tailed investigaiion  of  Lhc  synoptic  patterns  of  the  North  Pacific  Ocean, 
The  following  modifications  were  made  (Figs.   18,   19  and  20): 

1.  Vortex  B 

Vortex  B  is  poorly  formed  and  the  closed  center  has  dissipated. 
No  modifications  were  attempted  on  this  system. 

2.  Vortex  C 

Vortex  C  appears  diffuse.    The  500  mb  analysis  was  modified 
to  reflect  the  phase  and  amplitude  of  the  trough  lying  along  155W,  as 
indicated  by  the  open  area  to  the  west  and  enhanced  cloudiness  to  the 
east. 

3.  Vortex  D 

Vortex  D  has  continued  to  dissipate.    The  analysis  was   modi- 
fied only  to  accentuate  the  ridging  which  contributes  to  the  clear  area 

in  the  photo  noted  at  40N  between  vorticc;s  D  and  E. 
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4.  Vortex  E 

Vortex  E  is  moving  rapidly  to  the  east  northeast  across  the 
mid- Pacific  Ocean  at  45  kt.    The  cold  dry  tongue  has  changed  slightly 
but  the  system  has  not  intensified  appreciably.    It  has  deepened  to  992 
mb  at  sea  level  and,  based  on  the  ITOS  video  mosaic,  the  system's  posi- 
tion was  altered  slightly  westward  and  modified  to  obtain  the  closed  cir- 
culation and  a  500  mb  height  of  5160  m,  for  continuity. 

5.  Vortex   F 

Vortex  F  approaches  the  west  coast  of  North  America  as  a 
fully    mature  system  with  a  cloud  shield  that  covers  the  northwestern 
United  States  and  Canada.    The  analysis  was  modified  to  maintain  con- 
tinuity in  intensity  and  position  based  on  the  satellite  pictures. 

6 .  Vortex  G 

Vortex  G  is  dissipating  and  no  further  modifications  were 
attempted  on  this  system. 

7 .  Vortex  H 

Vortex  H  has  developed  rapidly  since  1200  GMT    11  March 
and  the  video  photo  clearly  displays  the  enhanced  comma  shaped  cloud 
(PVA  Max)  which  appears  just  ahead  of  the  base  of  a  short  wave  trough, 
indicated  by  the  cloud-free  subsiding  air  imbedded  in  the  long  wave 
ridge  as  it  interacts  with  the  frontal  band.    The  analysis  was  modified  to 
reflect  the  curvature  and  intensity  associated  with  this  secondary  frontal 
wave. 
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The  modification  effort  was  successful  with  this  short-wave  feature 
as  the  numerical  500  mb  MOD  ANAL  accepted  bogus  data  that  contoured 
the  trough  with  deepening  as  much  as  160  m  in  one  location. 

The  associated  ridge  to  the  west  was  built  into  the  trough  of  the 
UNMOD  ANAL  at  150W.    To  depict  the  increasing  baroclinicity  associated 
with  this  system,  the  gradient  was  tightened  in  the  southeast  quadrant  of 
vortex  H  as  far  east  as  125W. 
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FIG.     7 
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'ig.    8 


Fig.  6  500  mb  analyses 

0000  GMT  10  March  1971 
Values  are  5XX0  meters 


UNMOD  ANAL 

MOD  ANAL 

(no  modification  made) 


Fig„  7  Sea  level  analyses 

0000  GMT  10  March  1971 
Values  are  10XX  or  9XX  mb 


UNMOD  ANAL 
MOD  ANAL 


Fig.  8  ITOS-1  Video  Mosaic  for  approximately  0000  GMT  10  March 

1971  at   140W. 
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Fig.   11 


Fig.  9  500  mb  analyses 

1200  GMT  10  March  1971 
Values  are  5XX0  meters 


UNMOD  ANAL 
MOD  ANAL 


Fig.   10        Sea  level  analyses 

1200  GMT  10  March  1971 
Values  are  10XX  or  9XX  mb 


UNMOD  ANAL 
MOD  ANAL 


Fig.   11       Photo  of  ITOS-IR  Mosaic  for  approximately  1200  GMT  10  March 
1971  at   140W. 
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FIG.    13 
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Fig.    14 


Fig.   12        500  mb  analyses 

0000  GMT  11  March  1971 
Values  are  5XX0  meters 


UNMOD  ANAL 
MOD  ANAL 


Fig.  13        Sea  level  analyses 

0000  GMT  11  March  1971 
Values  are  10XX  or  9XX  mb 


UNMOD  ANAL 
MOD  ANAL 


Fig.  14       ITOS-1  Video  Mosaic  for  approximately  0000  GMT  11  March 
1971  at  140W. 
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FIG.    16 
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Fig.   17 


Fig.   15        500  mb  analyses 

1200  GMT  11  March  1971 
Values  are  5XX0  meters 


UNMOD  ANAL 
MOD  ANAL 


Fig.  16        Sea  level  analyses 

1200  GMT  11  March  1971 
Values  are  10XX  or  9XX  mb 


UNMOD  ANAL 
MOD  ANAL 


Fig.   17        Photo  of  ITOS-IR  Mosaic  for  approximately  1200  GMT  11  March 
1971  at  140W. 
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FIG.    19 
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Fig.   20 


Fig.   18        500  mb  analyses 

0000  GMT  12  March  1971 
Values  are  5XX0  meters 


UNMOD  ANAL 
MOD  ANAL 


Fig.   19        Sea  level  analyses 

0000  GMT  12  March  1971 
Values  are  10XX  or  9XX  mb 


UNMOD  ANAL 
MOD  ANAL 


Fig.  20       ITOS-1  Video  Mosaic  for  approximately  0000  GMT  12  March 
1971  at  HOW. 
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VI.      RESULTS 

To  facilitate  the  discussion,  the  experiments  in  modifying  the 
analyses  and  the  effect  on  the  attendant  prognoses  will  be  taken  up 
separately.    Comments  concerning  the  utilization  of  the  four  techniques 
that  resulted  in  the  quantification  of  the  satellite  data  will  be  included 
in  the  analysis  discussion. 

A.         ANALYSIS 

The  results  of  the  experimental  500  mb  and  sea  level  pressure  modi- 
fications of  the  analyses  (MOD  ANALS)  from  0000  GMT  10  March  through 
0000  GMT  12  March  1971.-  and  the  corresponding  analyses  from  conven- 
tional data  (UNMOD  ANALS)  are  displayed  opposite  the  applicable  satel- 
lite pictures  in  section  V.      It  must  be  pointed  out  that  the  desired  degree 
of  modification  (see  Appendix  C)  could  not  be  objectively  (numerically) 
reproduced  in  all  cases.    However,  only  the  objectively  modified  analyses 
are  presented  in  section  V. 

The  influence  of  the  satellite  data  on  the  analyses  is  readily  appar- 
ent in  each  of  the  figures.    As  an  example,  Fig.   10  shows  what  is  con- 
sidered to  be  the  most  pronounced  case  of  satellite  modification  made  to 
the  objective  upper-air  analyses  for  the  particular  period  of  time.   The 
modification  changes  the  analysis  of  zonal  flow  to  a  more  meridional  flow 
around  short-wave  troughs  and  ridges. 
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The  rapid  intensification  of  Vortex  F  between  10  and  11  March  1971 
is  an  outstanding  example  of  the  situation  encountered  by  the  analyst 
(human  or  machine)  working  in  the  data  sparse  area  of  the  North  Pacific 
Ocean.    Using  the  FWC  Pearl  Harbor  manually  analyzed  sea-level  pres- 
sure maps  purely  as  an  example,  vortex  F,  when  measured  from  the  erro- 
neously positioned  peak  of  the  frontal  wave  at  1200  GMT  10  March  to  its 
accurately  specified  position  12  hours  later  seemingly  moved  at  the  exces- 
sive speed  of  80  kts  across  much-used  shipping  lanes,  dropped  34  mb  in 
pressure,  and  increased  its  wind  from  10-15  kts  to  a  storm-force  intensity 
of  55  kt.    Real-time  availability  of  satellite  photos  or  a  nephanalysis  at 
0000  GMT  and  application  of  the  satellite  forecasting  rule  —  that  a  PVA 
maximum  approaching  a  convex  bulge  in  a  trailing  cold  front  may  result  in 
a  rapidly  deepening  surface  cyclone  —  would  have  permitted  a  re-analysis 
of  the  position  of  the  frontal  wave  and  its  potential  for  deepening.    The 
photo  of  excellent  ITOS-1   IR  nighttime  observations  at  1246  GMT  10  March 
would  have  confirmed  the  presence  of  a  deepening  system,  and  the  first 
daylight  ATS-1  and  ESSA-8  video  pictures  about  1800  GMT  10  March  would 
have  emphasized  the  need  for  a  storm  warning  superseding  previously 
issued  gale  warnings  or  heavy  weather  advisories. 

Comments  concerning  the  four  data  input  methods  of   1)  subjective 
techniques,    2)  vector  winds ,    3)  the  Nagle-Hayden  technique,  and 
4)   SIRS  soundings, which  were  used  in  the  modifications  follow. 
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1 .         Subjective  Techniques 

Subjective  methods  of  incorporating  satellite  data  entered 
the  modification  efforts  in  a  number  of  ways: 

a.  Revising  Conventional  Soundings 

A  revision  of  rejected  and  missing  portions  of  conven- 
tional rawinsonde  observations  through  comparison  with  satellite  photos 
and  SIRS  soundings  led  to  several  significant  improvements  of  the  con- 
ventional data.     For  example,  the  Ship  Papa  500  mb  wind  at  1200  GMT 
10  March  1971  failed  the  FNWC  Monterey  consistency  check  and  was 
rejected.      (See  1200  GMT  10  March  modification  discussion.)      It  was  a 
simple  matter  to  compare  the  cloud  patterns  in  the  vicinity  of  Ship  Papa 
with  the  winds  reported  on  the  sounding  above  and  below  the  500  mb 
level,  determine  the  probable  transmission  garble,  and  correct  it  with  a 
bogus  wind  report. 

b.  Determining  Intensification  of  Synoptic  Features 

Two  basic  steps  were  followed  in  determining  the  amount 
of  intensification  undergone  by  a  given  synoptic  feature  from  one  map  time 
to  the  next:     1)  close  surveillance,  and   2)  choosing  a  satellite  photogra- 
phy analog. 

Close  surveillance  was  afforded  by  the  outstanding 
collection  of  satellite  photographs  obtained  from  the  various  agencies 
noted  in  Table  1 .      The  time  lapse  in  photographic  coverage  was  never 
greater  than  12  hours.    The  photographs  of  nighttime  IR  values  were  parti- 
cularly valuable  since  these  mosaics  nearly  coincided  with  the  1200  GMT 
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synoptic  time.    Consequently,  during  the  time  frame  of  the  study  each 
synoptic  feature,  particularly  the  vortical  systems,  could  be  clearly 
followed. 

Analogs  were  selected  from  the  various  satellite  photog- 
raphy guidebooks    [1,  7,   12,    14,   15,   17  and  26]       which  contain  an  excel- 
lent assortment  of  time-lapse  photographs  relating  changes  in  cloud  patterns 
to  synoptic  map  changes  for  the  given  feature  (i.e.    wave  cyclone,  cold 
trough,  etc.). 

Using  the  analog  selected,  and  noting  the  changes  in 
the  individual  cloud  patterns,  a  qualitative  decision  was  made  concerning 
the  sign  of  the  intensity  change  since  the  last  map  time.    The  confidence 
level  of  the  decision  was -strengthened  by  "looking  ahead"  to  the  next  0600 
or  1800  GMT  sea-level  pressure  analysis  (see  section  IV  for  an  explana- 
tion of  this  "updating"  procedure).     Next,  tentative  500  mb  height  values 
were  assigned  to  the  feature  on  the  basis  of  the  12-hour  old  value,  the 
current  operational  analysis  value,  and  the  subjectively  indicated  change. 
This  tentative  value  was  then  compared  with  the  computation  from  the 
Nagle-Hayden  technique  and  the  SIRS  values  (if  available)  for  a  "best  fit" 
with  the  conservative  thickness  pattern.    This  final  value  was  then  placed 
into  the  modification. 

c.         Tightening  of  Gradient  in  Proximity  to  Jet  Streams 

The  500  mb  gradient  was  geostrophically  corrected  in 
the  vicinity  of  the  polar  jet  stream  where  clear-cut  detection  of  jet  stream 
cirrus  was  possible     [15  ]  .    This  method  would  have  been  more  useful  for 
modifying  the  300  or  200  mb  levels. 
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Several  excellent  satellite  photographs  of  subtropical 
jet  streams  were  noted  during  this  period,  but  no  attempt  was  made  to 
modify  the  wind  flow  in  the  lower  latitudes . 

d.         Modifying  Thickness  Contours  near  Fronts 

After  surface  fronts  were  re-positioned  on  the  basis  of 
satellite  interpretation,  the  1000/500  mb  thickness  patterns  were  adjusted 
to  conform  to  the  standard  model  pattern  of  the  National  Weather  Service 
[25  ]  .       After  revising  the  surface  analysis,  if  necessary,  the  500  mb 
field  was  retrieved  by  graphical  addition  of  the  surface  and  thickness 
fields.    This  procedure  proved  quite  simple  and  rapid  to  use  in  practice 
since  the  area  of  change  was  normally  small.     Furthermore,  when  used  in 
conjunction  with  the  Nagle-Hayden  SD  zero-vorticity  line  delineation, 
this  technique  invariably  resulted  in  a  more  meridional  trough  to  ridge 
flow  than  before  at  500  mb  over  the  associated  surface  cyclone  feature. 
2.         Vector  Winds 

The  ATS-1  wind  reports  recorded  at  NESS  by  computing  the 
vector  movement  of  cloud  patterns  recorded  on  film  loops  did  not  provide 
useful  inputs  during  this  particular  period.    The  migratory  pressure  patterns 
were  generally  north  of  4 ON,  whereas  vector  winds  were  seldom  computed 
north  of  30N.    These  computations  were  restricted  to  the  300  mb  level, 
though  NESS  does  derive  low-level  (850  mb)  winds  on  occasion.    Unques- 
tionably, low-level  wind  observations  during  the  development  of  vortex  F 
would  have  been  extremely  useful.    Though  the  authors  had  access  to  the 

o 

film  loops,  there  was  no  vector  equipment  on  hand. 


8  Vector  wind  console  and  satellite  picture  equipment  are  scheduled  for 
installation  at  the  FNWC  Monterey  by  January  1972. 
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3.  Nagle-Hayden  Technique 

The  Nagle-Hayden  technique  for  deriving  500  mb  height  values 
from  the  SR  500  mb  field  by  using  the  geometry  of  certain  types  of  observ- 
able  cloud  patterns  proved  extremely  workable  in  practice.    Once  the 
authors  became  accustomed  to  the  methods  of  delineating  the  SD  zero- 
line  and  spiral  cloud  amplitude  and  diameter,  the  corresponding  statisti- 
cal heights  for  the  North  Pacific  Ocean  could  be  constructed  in  less  than 
one  hour.    Presumably,  a  skilled  analyst  could  perform  this  task  in  less 
than  30  minutes.    With  access  to  a  CRT  (cathode  ray  tube)  console  at  a 
major  computer  center,  the  entire  process  is  amenable  to  computerization, 
with  the  values  being  derived  and  incorporated  automatically  in  a  matter 
of  minutes . 

The  zero-line  technique  was  particularly  useful  in  interpreting 
short-wavelength  troughs  and  ridges.    At  each  map  time,  this  method  was 
used  to  derive  500  mb  heights  along  troughs  and  ridges  across  the  entire 
ocean.    Although  some  of  the  values  were  then  changed  subjectively, 
most  of  them  could  have  been  incorporated  without  change  into  the  MOD 
ANAL  if  the  FNWC  restriction  of  50  bogus  points  could  have  been  expanded. 

4.  SIRS  Data 

The  FNWC  program  restriction  of  fifty  bogus  points  at  the  500 
mb  level  prohibited  direct  input  of  all  the  available  SIRS  soundings .    The 
SIRS  values  for  500  mb  were  considered  as  data  during  the  manual  MOD 
ANAL  effort  and  therefore  did  influence  the  bogus  values  used  to  achieve 
the  desired  objective  MOD  ANAL.    Available  SIRS  coverage  for  the  North 
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Pacific  Ocean  was  poor  on  the  1200  GMT  10  March  map  (two  reports), 
fair  at  0000  GMT  11  March  (11  reports)  and  1200  GMT  11  March  (10  re- 
ports), and  good  at  0000  GMT  10  March  (19  reports)  and  0000  GMT  12 
March  (21  reports).     Of  these  soundings,  many  were  located  in  the  sub- 
tropical regions,  and  not  used  in  the  modification.    The  remaining  sound- 
ings were  useful  in  the  re-positioning  of  troughs  and  ridges. 

To  cite  an  example,  on  a  north-south  oriented  trough  line  along 
155W  at  0000  GMT  10  March,  three  SIRS  heights  of  5190  m,   5370  m,  and 
5680  m  were  recorded.    The  corresponding  FNWC  operational  analysis 
values  were  5310  m,   5520  m,  and  5760  m.    The  input  of  the  SIRS  values 
would  have  improved  the  FNWC  operational  analysis  of  this  short-wave 
trough.    At  0000  GMT  12  March,  in  the  vicinity  of  a  deepening  trough 
near  135W,  two  SIRS  values  of  5490  m  and  5480  m  compared  with  5640  m 
and  5620  m,  respectively,  for  the  FNWC  analysis. 

In  summary,  the  scanty  amount  of  applicable  SIRS  500  mb  height 
values  was    used  to  great  advantage  in  the  re-positioning  of  short-wave 
troughs  and  ridges.    The  hemispheric  coverage  of  SIRS  500  mb  heights  for 
this  period  is  discussed  in  Appendix  B. 

B.         PROGNOSES 

The  FNWC  primitive  equation   (PE)  prediction  model  was  used  to 
generate  12,  24,  36,  and  48  hour  prognostic  fields  from  the  MOD  and 
UNMOD  ANALS  of  1200  GMT  10  March  and  0000  GMT  11  March  1971. 
The  sea   level  and   500  mb  MOD  and  UNMOD  PROGS   for  each  of  the 
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four  forecast  intervals  are  displayed  in  Figs.  21  through  36,  with  MOD 
(solid  lines)  and  UNMOD  (dashed  lines)  versions  superimposed. 
1 .         Droad  Aspects 

Salient  features  of  the  prognoses  are  summarized  below. 

1)  The  MOD  and  UNMOD  prognoses  appear  to  converge  to 

a  similar  pattern  with  increasing  forecast  interval;  this  fact  becomes  quite 
evident  by  the  48-hour  prognosis.    That  is,  the  small-scale  features  ap- 
pear to  be  smoothed  out  with  time  such  that  very  similar  large-scale 
patterns  remain  in  both  MOD  and  UNMOD  versions  by  48  hours. 

2)  Evidently,  the  FNWC    PE  model  is  very  stable  in  the 
matter  of  handling  bogused  short  wave  features.     No  large  errors  were 
generated  in  any  of  the  prognoses.    Dispersion  of  energy  from  the  new 
short  wave  features  across  the  Northern  Hemisphere  is  minimal.    In  fact, 
the  comparable    48  hour  prognoses  were  virtually  identical  in  the  regions 
not  affected  by  the  apparent  propagation  of  the  modified  synoptic  features. 
Fig.  37  illustrates  this  stability  by  showing  a  field  of  differences  between 

the  error  fields  of  the  48  hour  500  mb  MOD  and  UNMOD  prognoses.    The 

9 
error  fields  were  generated  by  verifying  with  the  500  mb  STACK  ANAL  (SA). 

That  is, 

Diff .  Field  =    [  (MOD  PROG  -  SA)  -  (UNMOD  PROG  -  SA)  ] ,  (6) 

which  is  the  same  as  taking  the  difference  between  the  two  prognoses. 
Small  differences  in  this  field  outside  the  area  of  interest  indicate  negli- 
gible side  effects  of  using  bogus  data  in  the  FNWC  PE  model. 


9STACK  ANAL  =  operational  FNWC  update  analysis. 


69 


2 .         The  Prognosis  of  Individual  Synoptic  Features 

a.         Discussion  of  1200  GMT  10  March  1971  Prognoses 
(Figs.  21  through  28). 

(1)  Vortex  B.    The  MOD  PROG  moved  vortex  B  very 
slowly  to  the  northwest,  and  filled  it  4  mb  to  a  976  mb  sea-level  value 
in  48  hours.    At  500  mb  there  was  only  slight  movement  northward.    The 
center  dropped  30  m  to  5050  m  in  the  first  12  hours  and  remained  at  that 
height  for  the  rest  of  the  48  hours. 

In  the  UNMOD  PROG  the  movement  was  slightly 
faster  but  the  surface  low  remained  in  the  Gulf  of  Alaska  and  deepened 
the  center  from  982  mb  to  974  mb  in  48  hours.    A  secondary  cyclone  is  gen- 
erated at  the  surface  within  12  hours  to  the  southeast  of  vortex  B  by  both 
the  MOD  and  UNMOD  prognoses.    The  24  hour  MOD  PROG  moved  this  low 
inland  and  filled  it.    The  UNMOD  PROG  retained  this  low  and  kept  it 
quasi-stationary  west  of  Vancouver  Island,  B.  C.  ,  out  to  48  hours.    Its 
central  pressure  dropped  from  9  89  mb  to  978  mb.      This  UNMOD  stationary 
feature  coincided  with  the  actual  position  of  the  fast  moving  vortex  F  at 
1200  GMT  12  March  and  resulted  in  the  4  8  hour  UNMOD  PROG  verifying 
well  in  this  area.    However,  its  quasi-stationary  position  during  this 
prognosis  is  not  verified  by  the  spiral  cloud  patterns  for  this  period  of 
time . 

(2)  Vortex  C.    The  cloud  vortex  is  represented  on  the 
500  mb  MOD  ANAL  as  a  cyclone;  its  corresponding  surface  low  repositioned 
westward  300  nmi  for  proper  vertical  alignment.    The  MOD  PROG  moved 
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this  vortex  eastward  at  30  kt  for  the  first  12  hours,  and  thereafter  it 
averaged  20  kt  out  to  48  hours.    Slight  filling  occurred  in  the  first  12  hours 
but  the  closed  system  was  retained  at  the  surface  out  to  48  hours,  having 
filled  3  mb  at  sea  level  and  130  m  at  500  mb,  where  the  trough  is  in  evi- 
dence at  155W. 

The  UNMOD  PROG  moved  this  system  eastward 
at  48  kt  and  filled  it  to  where  only  troughs  remained  at  both  sea  level  and 
500  mb  by  36  hours,  followed  then  by  further  damping  of  amplitude  to  48 
hours . 

(3)  Vortex  D.    This  small,  cold-core  system  is 
moved  slowly  eastward  and  filled  by  both  MOD  and  UNMOD  PROGS. 

(4)  Vortex  E.    This  system  is  modified  in  the  analysis 
at  500  mb  and  given  a  better  trough-ridge  configuration  with  a  60  m  height 
change.    MOD  and  UNMOD  PROGS  handled  the  movement  of  the  sea-level 
system  similarly  with  a  northeastward  movement  at  13  kt.    However,  at 
the  500  mb  level,  only  the  MOD  PROG  retained  this  short-wave  feature 
out  to  36  hours.    The  UNMOD  PROG  version  assimilated  this  feature  after 
12  hours . 

(5)  Vortex  F.    This  closed  circulation  was  placed  in- 
to an  existing  sea  level  trough  on  the  MOD  ANAL  but  was  lost  on  the  12 
hour  MOD  PROG.     However,  the  trough  was  retained  in  the  MOD  PROG 
but  not  in  the  UNMOD  PROG.      On  the  MOD  PROG  the  trough  moves  to  the 
northeast  and  deepens  into  a  closed    circulation  in  48  hours.    Though  all 
the  subsequent  cloud  patterns  indicate  this  moving  trough  was  continuously 
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a  closed  circulation,  the  movement  in  the  MOD  PROG  appears  good  and 
verifies  well  in  48  hours.    As  mentioned  earlier,  the  UNMOD  PROG 
missed  this  moving  system,  but  by  coincidence  the    quasi  stationary  sea- 
level  low  in  the  UNMOD  PROG  verifies  very  well  on  the  48  hour  prognosis 
(but  only  on  this  prognosis)  relative  to  the  observed  cloud  patterns. 

(6)        Vortex  G.    This  system  was  handled  identically 
and  fairly  well  by  both  MOD  and  UNMOD  PROGS.      It  was  moved  too 
slowly  at  500  mb  causing  the  cyclone  to  develop  some  vertical  tilt. 

b.         Discussion  of  0000  GMT  11  March  1971  Prognoses 
(Figs.   29  through  36). 

(1)  Vortex  B.      The  sea  level  MOD  and  UNMOD  PROGS 
moved  vortex  B  toward  the  north-northwest  (NNW)  with  the  UNMOD  PROG 
fluctuating  around  975  mb  and  the  MOD  PROG  filling  the  center  by  12  mb 
to  983  mb.    The  centers  coincide      in  the  northern  Gulf  of  Alaska  at  48 
hours.    As  cyclone  B  moves  towards  the  NNW,  the  troughs  to  the  south- 
east of  vortex  B  remained  stationary  on  both  the  sea  level  and  500  mb 
UNMOD  PROGS. 

(2)  Vortex  C.    The  MOD  PROG  sea  level/500  mb, 
982  mb/4970  m  low  moves  east-southeast  (ESE)  at  30  kt  and  fills  to  987 
mb/5040  m  in  12  hours  and  slows  to  20  kt  while  filling  to  991  mb/5220  m 
in  48  hours.    To  the  northeast  at  500  mb,  a  short-wave  ridge  appears  in 
the  long-wave  trough  in  24  hours,  but  is  absorbed  in  the  large-scale  flow 
by  48  hours. 
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The  corresponding  trough  in  the  UNMOD  PROG 
fills  rapidly  and  is  completely  absorbed  by  the  long-wave  flow  in  24  hours. 

(3)  Vortex  D.    This  small  system  is  handled  well  by 
both  prognoses  as  it  moved  slowly  eastward  and  dissipated. 

(4)  Vortex  E.    The  MOD  PROG  moves  this  system  at 
20  kt  toward  the  east-northeast  and  deepens  it  from  1010  mb  to  990  mb  at 
sea  level  in  48  hours.    The  associated  500  mb  trough  is  continually 
broadened  and  sustained  while  the  previously  enhanced  ridge  to  the  east 
is  completely  smoothed  into  the  long-wave  pattern. 

The  UNMOD  sea  level  trough  moved  eastward  at 
29  kt  and  the  circulation  intensified,  but  the  center  only  deepened  slightly. 
There  is  no  500  mb  trough  development  associated  with  the  sea  level  trough 
for  this  UNMOD  PROG. 

(5)  Vortex  F.      This  system  existed  only  in  the  MOD 
ANAL  and  PROG.     Considerable  filling  occurred  in  the  first  12  hours  with 
the  sea-level/500  mb  centers  changing  from  975  mb/5090  m  to  982  mb/ 
5160  m  and  retained  only  a  weak  trough.     By  24  hours  the  sea  level  trough 
progressed  northeast  and  blended  with  the  stationary  trough  southeast  of 
vortex  B.    Similarly  at  500  mb,  vortex  F  was  no  longer  present  as  a  dis- 
tinctive feature,  having  merged  into  the  long-wave  trough  dominating  the 

region. 

(6)  Vortex  G.    This  system  was  not  modified  and  both 
prognoses  handled  it  equally  well,    but  excessive  tilting  of  this  system 
occurred  with  increasing  forecast  interval. 
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FIG.  21  500  MS  PROGNOSES 
1200  GMT  10  MARCH  1971 
,  VALUES  ARE  5XX0  METERS 


--UNMOO  PROG 
—  MOD  PRQG 
2  Hi?  PROG 


FIG  22  SEA  LEVEL  PROGNOSES 
1200  GMT  10  MARCH  1971 
k     VALUES  ARE  10XX  OR  9XX  MB 


---UNMOD  PROG 
7— MOD  PROG 
2  HR  PROG 
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FIG.  23  500  M3  PROGNOSES 
1200  GMT  10  MARCH  1971 
*  VALUES  ARE  5XX0  METERS 


---UNMOD  PROG 
— MOD  PROG 
24  HR  PROG 


FIG  24  SEA  LEVEL  PROGNOSES 
1200  GMT  10  MARCH  197  1 
.  VALUES  ARE  10XX  OR  9XX  MB 


---UNMOD  PROG 
— MOD  PROG 
24  HR  PROG 
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FIG. 26  SEA  LEVEL  PROGNOSES 
1200  GMT  10  MARCH  1971 
*  VALUES  ARE  10XX  OR  SXX  MB 


---UNMOO  PROG 
— MOD  PROG 
36  HR  PROG 
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FIG. 27  500  MB  PROGNOSES 
1200  GMT  10  MARCH  1971 
*  VALUES  ARE  5XX0  METE.7S 


---UNMOD  PROG 
— MOD  PROG 
48  HR  PROG 


iJSn2?u?EinLEVEL  PROGNOSES 
J?00  GMT  10  MARCH  1971 

VALUES  ARE  10XX  OR  9XX  M3 


---UNMOD  PROG 
— MOD  PROG 
48  HR  PROG 
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FIG. 29  500  MB  PROGNOSES 
.  0000  GMT  11  MARCH  1971 
>  VALUES  ARE  5XX0  METERS 


---UNMOD  PROG 
— MOD  PROG 
12  HI?  PROG 


.FIG  30  SEA  LEVEL  PROGNOSES 
0000  GMT  11  MARCH  1971 
VALUES  ARE  10XX  OR  SXX  MB 


— MOD  PROG 
12  HR  PROG 
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FIG. 31  500  M3  PROGNOSES 
0000  GMT  11  MARCH  1971 
a  , VALUES  ARE  5XX0  METERS 


---UNMOD  PROG 
— MOD  PROG 
24  HR  PROG 


FIG. 32  SEA  LEVEL  PROGNOSES 
0000  GMT  11  MARCH  1971 
*'  VALUES  ARE  10XX  OR  9XX  M3 


---UNMOD  PROG 
—  MOD  PROG 
24  HR  PROG 
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FIG. 33  500  MS  PROGNOSES 
0000  GMT  11  MARCH  1971 
VALUES  ARE  5XX0  METERS 


---UNMOD  PROG 
— MOD  PROG 
36  HR  PROG 


FIG. 34  SEA  LEVEL  PROGNOSES 
0000  GMT  11  MARCH  1971 
VALUES  ARE  10XX  OR  9XX  MB 


---UNMOD  PROG 
— MOD  PROG 
36  HR  PROG 


80 


FIG. 35  500  MB  PROGNOSES 
IGOOO  GMT  1  I  MARCH  1 971 
'VALUES  *nE  5XX0  METERS 


---UNMOD  PROG 
— MOD  PROG 
•48  H3  PPOG 


FIG. 3 6  SEA  LEVEL  PROGNOSES 
0000  GMT  11  MARCH  197 
VALUES  ARE  10XX  OR  9XX  MB 


---UNMOD  PROG 
— MOD  PROG 
48  HR  PROG 
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FIG.    37   500. MB. 


PROJECTION     POLAR  STEREOGRAPHIC-TRUE  AT  60     NORTH   LATITUDE 
SCALE     1.60.000.000 


FLEET  NUMERICAL  WEATHER  CENTRAL 
MONTEREV,  CALIFORNIA 


CHART  No  69- 1 


Fig.  37      500  rnb  48  hour  Prognoses  Differential  Field  obtained  from  the 
MOD  and  UNMOD  PROG  verification  against  the  STACK  ANAL. 
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C.         VERIFICATIONS 

The  problem  of  selecting  a  suitable  verification  procedure  for  satel- 
lite-derived modifications  has  been  stressed  by  many  fellow  workers 
[2,   4,   8]    .      This  study  emphasizes  the  use  of  satellite  data  for  modifi- 
cations to  synoptic  analysis  over  the  sparse-data  oceanic  regions  where 
cloud  patterns  are  essentially  void  of  complicating  topographic  effects. 

The  first  part  of  the  problem  involves  obtaining  an  accurate  analysis 
of  the  oceanic  atmosphere  taken  only  from  a  broad  base  of  conventional 
data  over  the  North  Pacific  Ocean,  against  which  the  MOD  PROGS  and 
the  UNMOD  PROGS  could  be  verified.    The  authors  resolved  this  part  of 
the  problem  by  using  the  FNWC  update  analyses,  henceforth  referred  to 
as  the  STACK  ANALS  which  are  largely  free  of  satellite  input. 

The  second  part  of  the  problem  is  the  matter  of  obtaining  the  most 
accurate  analysis  using  the  maximum  of  conventional  and  satellite  data 
for  the  basis  of  the  analyses,  against  which  only  the  MOD  PROGS  would 
be  verified.    This  part  of  the  problem  is  resolved  by  using  the  authors' 
MOD  ANALS. 

1 .         Description  of  Prognostic  Verification  Schemes  Used 

The  initial  evaluation  of  the  prognoses  is  the  difference  fields 
between  the  MOD  PROG  and  verifying  STACK  ANAL  and  the  UNMOD  PROG 
and  the  same  STACK  ANAL,  as  MOD  PROG-STACK  ANAL  and  UNMOD  PROG- 
STACK  ANAL  (Figs.  38  through  53).    With  this  comparison,  the  effects  of 
initialization  from  both  the  modified  and  unmodified  fields  on  the  resulting 
prognoses  will  be  measured  against  the  best  available  verifying  FNWC 

analyses. 
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The  second  evaluation  of  prognostic  results  considers  a  comparison 
of  the  following  difference  fields,   MOD  PROG- STACK  ANAL  and  UNMOD 
PROG -MOD  ANAL  (Figs.   54  through  63). 

Thus,  the  first  evaluation  assumes  the  STACK  ANALS  are  the  best 
available  for  verification  while  the  second  evaluation  selects  the  best 
available  corresponding  analyses  for  verification  of  each  prognosis. 

2 .         Summary  of  Difference  Field  Evaluations 

A  study  was  performed  on  all  difference  fields  generated.     The 
patterns  were  evaluated  and  are  discussed  in  three  contexts.    Absolute 
differences  equivalent  to  or  exceeding  4  mb  at  sea  level  and  60  m  at 
500  mb  are  considered  significant  for  this  evaluation.    This  criteria  im- 
poses a  closer  tolerance  at  sea  level  than  at  500  mb. 
a.         General  Overall  Evaluation 

The  following  comments  apply  to  MOD  and  UNMOD 
PROGS  in  relation  to  both  STACK  and  MOD  ANALS. 

The  relative  locations  of  major  (two  or  more  closed 
contours)  difference  centers  dramatically  and  consistently  reflect    a  high 
correlation  with  the  position  of  analyzed  low  centers,  troughs  or  ridges. 
(There  were  no  major  high  centers  in  the  study  area.)    In  particular,  the 
locations  of  the  major  negative-difference  centers  correlate  with  the 
analyzed  positions  of  ridges  while  the  locations  of  the  major  positive- 
difference  centers  correlate  with  the  analyzed  positions  of  low  centers 
and  troughs.    This  marked  correlation  strongly  suggests  that  the  short- 
wave-amplitudes in  the  prognoses  have  been  damped.    Wave-amplitude 
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damping  of  the  short-wave  features  dominates  the  contributions  to  the 
difference  patterns  and,  therefore,  masks  some  phase  difference  contri- 
butions for  these  waves.    However,  there  are  detectable  phase  problems 
in  the  short-wave  pattern  which  result  from  the  prognostic  model  moving 
the  short-wave  trough-ridge  systems  too  slowly  compared  to  the  verifying 
analyses.    This  contribution  appears  to  increase  with  increasing  forecast 
interval  but  remains   a    relatively  minor    contribution  to  the  total  differ- 


ence field . 


In  addition,  a  qualitative  evaluation  shows  that  the 


magnitudes  and  areal  coverage  of  major  difference  features  are  similar  in 
the  UNMOD  and  MOD  PROG  verifications. 

b.         Compiled  Difference  Field  Results 

Every  difference  field  for  which  both  MOD  and  STACK 
ANALS  were  available  is  compared  and  quasi-objectively  evaluated  on  the 
basis  of  apparent  forecast  accuracy.    The  results,  listing  the  prognoses 


that  verify  best,  are  shown  in  Table  II. 

PROG  from 

1200  GMT    10  Mar 


PROG  from 

0000  GMT  11  Mar 


PROG  II 

verifie: 

4TERVAL 
D  AGAINST 

12 

24 

36 

12 

24 

STACK 
ANAL 

500  mb 

UNMOD 

MOD 

UNMOD 

UNMOD 

UNMOD 

Sea  Level 

UNMOD 

MOD 

EQUAL 

MOD 

MOD 

BEST 
CORRES 
ANAL 

500  mb 

MOD 

UNMOD 

UNMOD 

UNMOD 

UNMOD 

Sea  Level 

UNMOD 

MOD 

EQUAL 

MOD 

MOD 

Table  II.      The  results  of  a  quasi-objective  evaluation  of  all  difference 
fields  and  selecting  the  best  verification. 


Interpreting  the  results   of  this   evaluation  indicates 
that  the  modifications   improved  the  verification  of  the  majority  of  sea 
level  prognoses   but  did  not  improve  the  verification  of  the  majority  of 
the   500  mb  prognoses   against  the  STACK  ANALS . 

c.         Salient  Features 

(1)  Consistency  in  the  Difference  Patterns.    Relative 
to  the  verifying  analyses,  there  is  a  consistent  tendency  for  the  PE  model 
to  underforecast  the  amplitude  of  ridges  and  troughs ,  especially  the  short- 
wave features.    Since  this  type  of  error  becomes  evident  in  the  first  12 
hours  of  prognosis,  the  major  part  of  this  problem  may  be  related  to  the 
initialization  of  the  PE  model  rather  than  its  ability  to  handle  short-wave 
fsature  s 

(2)  Inconsistencies  in  the  Difference  Patterns. 
Possibly  relating  to  initialization  problems ,  the  amount  of  difference  occur- 
ring during  the  first  twelve  hours  of  prognosis  was  greater  than  the  differ- 
ence change  during  any  other  twelve  hour  interval.    However,  there  is 
little  consistency  in  the  amount  or  sense  of  the  differences  generated  during 
any  twelve  hour  period.     Likewise,  there  is  no  consistent  pattern  of  gener- 
ation of  significant  differences  for  a  particular  forecast  interval  (as  12, 

24,  36,  or  48  hours). 

(3)  Effects  of  Small-Scale  Features.    The  difference 
fields  reinforce  the  point  noted  in  section  VI.  B,  namely  that  many  small- 
scale  features  were  highly  smoothed  in  12  hours  and  completely  obliterated 
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by  24  hours.    This  lack  of  ability  to  retain  and  develop  the  small-scale 
ridges  and  troughs  accounts  for  a  significant  portion  of  each  of  the  MOD 
PROG  -  MOD  ANAL  difference  fields  and  neutralizes  much  of  the  effort  made 
to  modify  the  UNMOD  ANALS.     The  small-scale  features  were  the  major 
modifications  made  to  the  analyses,  especially  at  500  mb.    However,  it 
should  be  emphasized  that  these  small-scale  features  must  still  be  re- 
garded as  synoptic  scale  in  character  and,  therefore,  desirable  to  include 
in  analyses  or  prognoses. 

(4)        Effects  of  Using  Bogused  Reports.    There  are  no 
indications  that  the  use  of  bogus  values  to  modify  the  analyses  generates 
spurious  errors  in  the  prognoses,  either  in  the  general  area  of  the  initial 
modifications   ,  or  over  the  entire  hemisphere,  as  seen  in  Fig.  37. 
3 .         Comments  Relating  to  the  Difference  fields 

Many  aspects  of  a  PE  model,  e.g.  ,  performance  characteristics, 
formulations,  and  the  assumptions  made,  are  understood  most  clearly  by 
those  who  designed  the  model.     However,  in  some  instances  the  model's 
behavior  has  not  been  fully  documented.    It  is  undoubtedly  true  that  the 
forecaster  in  the  field  has  the  least  insight  into  the  complex  non-linear 
interactions  which  take  place  in  a  primitive  equation  forecast    [11  ]. 

In  view  of  the  above,  this  section  is  directed  toward  an  in- 
creased understanding  of  the  utility  of  satellite  observations  in  relation  to 
its  potential  in  increasing  the  accuracy  of  forecasts  from  PE  models. 
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(1)  Causes  Related  to  the  PE  Model  and  Initialization 
Process.     Similar  systematic  error-pattern  features  from  NMC  PE  prognos- 
tics were  reported  by  Fawcett  in  1969  [3]  .    He  concluded  that  a  consistent 
pattern  of  positive  errors  in  monthly  mean  troughs  and  negative  errors  in 
monthly  mean  ridges  does  occur  on  the  average.    Further,  this  type  of 
error  pattern  shows  up  in  the  daily  prognostic  charts  as  a  lack  of  ampli- 
tude in  short-wave  troughs  and  ridges  and  an  underforecast  of  the  strength 
of  the  polar  jet  core.     More  recently  a  National  Weather  Service,  Western 
Region,  study  pointed  out  the  existence  of  "slowness  of  short-wave  trough 
motion  found  in  PE  forecasts , "    till    . 

The  existence  of  both  the  damping  and  slowness 
of  short-wave  features  has  been  noted  in  the  500  mb  prognoses  generated 
for  this  study.    The  significance  of  these  model  characteristics  with  re- 
gard to  the  analysis  modifications  applied  to  short-wave  features  cannot 
be  explained  in  this  research  because  the  initialized  fields  were  not  re- 
tained.   However,  it  is  felt  that  the  findings  in  this  study  corroborate  the 
findings  of  the  other  recent  aforementioned  studies  of  the  PE  model  [3,  11]  . 
It  is  stated  in   [11]    that,   "In  going  from  the  analysis  to  the  initialization, 
30  percent  of  the  energy  may  be  lost." 

(2)  Causes  Related  to  the  Analysis  and  Bogus  Tech- 
niques.     The  numerical  damping  and  retarding  of  the  short-wave  features 
prompted  further  study  of  the  manually  produced  MOD  ANALS  and  the  final 
numerically  generated  MOD  ANALS.    The  numerical  MOD  ANALS,  which 
were  accepted  by  the  authors ,  were  achieved  only  after  numerous  attempts 


at  selecting  the  best  bogus  values  and  locations  which  would  effectively 
result     in  the  numerical  MOD  ANAL  being  in  reasonable  agreement  with  the 
manual  MOD  ANAL.    It  is  noted  that  some  inconsistencies  do  exist  between 
the  positions  of  the  sea-level  and  500  mb  ridge  lines  in  the  final  numeri- 
cal MOD  ANALS.    This  condition  resulted  from  concentrating  primary  atten- 
tion on  selecting  bogus  values  to  properly  locate  troughs.    It  is  difficult 
to  completely  assess  the  total  effect  of  these  inconsistencies  on  the 
resulting  prognoses  of  short-wave  features. 

Another  unknown  factor  is  the  effect  of  not  explicit- 
ly modifying  the  soundings  above  500  mb.    In  effect,  the  FNWC  analysis 
program  extrapolates  both  the  500  mb  conventional  and  bogus  data  to  300 
mb  to  be  considered  along  with  all  other  data  for  analysis  at  the  latter 
level.    From  the  300  mb  analysis,  grid  point  values  are  then  extrapolated 
back  to  500  mb  to  aid  in  the  final  500  mb  analysis.    Thus,  the  300  mb 
analysis  may  affect  the  500  mb  analysis  in  such  a  way  as  to  partially  null- 
ify the  effect  of  the  500  mb  bogus  data.    In  summary,  the  effect  of  500  mb 
bogus  data  cannot  be  readily  determined  beforehand. 

The  lack  of  consistent  stratospheric  data,  con- 
ventional and  bogus ,  will  force  the  PE  model  to  climatological  values  which 
imply  stability  and  a  possible  damping  effect  on  the  desired  unstable  sys- 
tems which  could  assist  the  development  in  the  PE  prognosis. 
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FIG.  38  500  MB  VERIFICATIONS 

1200  GMT  10  MA3CH  1971  PROGNOSES 

12  HR  PROG  -  STACK  ANAL  IN  METERS.  VT  00  00  GMT  II  M4R 


iFIG  39  SEA  LEVEL  VERIFICATIONS        ---UMMOD 

1200  GMT  10  MARCH  1971  PROGiJOSES       Mm 

,  12  HR  PROG  -  STACK  ML  IN  Nl  VT  00  00  GMpHMAR 
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FIG.  40  500  M3  VERIFICATIONS  ---UNMGD 

1200  GMT  10  MA?CH  1971  PROGNOSES       — MOO 
[24  HR  PROG  -  STACK  ANAL  IN  METERS.  VT  1200  GMT  1 1  MAR 


FIG.  41  SEA  LEVEL  VERIFICATIONS 

1200  GMT  10  MARCH  1971  PROGNOSES       — MOD 

24  HR  PROG  -  STACK  ANAL  IN  M3-  VT  1200  GMT  11  MAR 
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FIG.  42    500  MB  VERIFICATIONS  ---UNMOD 

1200  GMT    10  MARCH   1971   PROGNOSES  —MOD 

36  HR  PROG   -   STACK   ANAL    IN  METERS.  VT  0000  GMT   T2  MAR 
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[FIG.    43  SEA  LEVEL  VERIFICATIONS  ---UNMOD  j 

1200  GMT    10  MARCH   1971   PROGNOSES  — MOD 

28  H3  PROG   -   STACK  ANAL    IN  M3  VT  0000  GMT    12  MAR 
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FIG.  4^  500  M3  VERIFICATIONS  ---UNMOD 

V2  00GMT  10  MARCH  1971  PROGNOSES       — MOO 
*  ;3  HR  PROG  -  STACK  ANAL  IN  METERS.  VT  1200  GMT  12  MAR 


FIG.  45  SEA  LEVEL  VERIFICATIONS 

12  00  GMT  10  MARCH  1S71  PROGNOSES       — 

■46  HR  PROG  -  STACK  ANAL  IN  M3.  VT  -I200  GMT  1 


UNMOD 
MOD 
2  MAR 
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FIG.  46  500  M3  VERIFICATIONS 
>00  00  GMT  11  MARCH  1971  PROGNOSES       — MOD 

V  12  HR  PROG  -  STACK  ANAL  IN  METERS.  VT  1200  GMT  11  MAS 


FIG.  47  SEA  LEVEL  VERIFICATIONS        ---IO0D 
O000  GMT  11  MARCH  1S71  PRGGilCSES       — MOD 
12  HR  PROG  -  STACK  ANAL  IN  M3.  VT  1200  GMT  11  MAR 


94 


FIG.    49  SEA  LEVEL   VERIFICATIONS 
;0000  GMT    11  MARCH   1971   PROGNOSES 
■   .24  HR  PROG  -  STACK  ANAL    111  M3.  VT  OOOO  GMT  12    MAR  - 
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FIG.  50  500  M3  VERIFICATIONS  ---UNMGD 

0000  GMT  11  MARCH  1971  PROGNOSES       — MCD 
v  36  HR  P.^OG  -  STACK  ANAL  IN  METERS.  VT  12  00  GMT  12  MAS 


FIG.  51  SEA  LEVEL  VERIFICATIONS 

0000  GMT  11  MARCH  1971  P^CGIICSES 

36  HR  PROG  -  STACK  ANAL  IN  M3.  VT  12  00  GMT  12  MAR 
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FIG.  52  500  M3  VER IF  I  CAT  I  CMS 
00  00  GMT  11  MARCH  1971  PROGNOSES 


a  48  HR  PROG  -  STACK  ANAL  IN  METERS. VT  0000  GMT  13  MAR 

I 


FIG.  53  SEA  LEVEL  VERIFICATIONS        — UNMOD 
00  00  GMT  II  MARCH  1971  PROGNOSES       —MOD 
A  48  HR  PROG  -  STACK  ANAL  IN  M3.  VT  0000  GMT  13  MAR 
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FIG.  54  500  MB  VERIFICATIONS 

1200  GMT  10  MARCH  1S71  PROGNOSES       - 

12  HR  PROG  -  BEST  C03RES  ANAL  IN  METERS.  VT  00  00  GMT  11  MAR 


FIG.  55  SEA  LEVEL  VERIFICATIONS        — UNMOD-STACK 

12  00  GMT  10  MARCH  1971  PROGNOSES       — MCD  -  MOD 
12  HR  PROG  -  DEST  CORRES  ANAL  IN  M3.  VT  00  00  GMT  11  MAR 
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FIG.  56  500  M3  VERIFICATIONS 

12  00  GMT  10  MARCH  1971  PROGNOSES 

OG  -  BEST  CORRES  ANAL  IM  METERS.  VT  12)0  GMT  11  MAR 


;FIG.  57  SEA  LEVEL  VERIFICATIONS 

'1200  GMT  10  MARCH  1971  PROGNOSES       - 

,£4  HR  PROG  -  DEST  CORRES  ANAL  IN  M3.  VT  12  00  GMT  11  MAR 
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FIG.  58  500  MB  VERIFICATIONS  ---UKMOD-STACK 

12  00  GMT  10  MARCH  1S71  P?CGi;OSES       — MOD  -  MOD 

36  HR  PROG  -  BEST  CORRES  ANAL  IN  METERS.  VT  0000  GMT  12  MAR 


FIG  59  SEA  LEVEL  VERIFICATIONS        ---UNMOD-STACK 
12  00  GMT  10  MARCH  1971  PROGNOSES       —MOD  -  MOD 
;  36  KR  PRCG  -  GEST  CORRES  ANAL  IN  MD.VTOO  00  GMT  12  MAR 
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FIG.  60  500  MS  VERIFICATION'S 

0000  GMT  11  MARCH  1S71  PROGNOSES 

12  HR  PROG  -  BEST  CORRES  ANAL  IN  METERS.  VT  I2  00  GMT  1 1  MAR 
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FIG.  61  SEA  LEVEL  VERIFICATION'S        — UNMOD-STACK 

OQOO  GMT  11  MARC!!  1971  PROGNOSES 

7?  HR  PROG  -  DEST  CORRES  ANAL  IN  MB.  VT  12  00  GMT  11  MAR 
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.   FIG.  62  500  MB  VERIFICATION'S  ---UNMOD-STACK 

0000  GMT  11  MARCH  1S71  PROGNOSES       — MCD  -  MOD 
24  HR  PROG  -  BEST  CORSES  ANAL  IN  METERS.  VT  00 00  GMT  12  MAR 


FIG.  63  SEA  LEVEL  VERIFICATION'S        ---UNMOD-STACK 
0000  GMT  11  MARCH  1971  PROGNOSES       — MCD  -  MCD 
k  24  HR  PROG  -  BEST  CORRES  ANAL  IN  MD.  VTOOOO  GMT  12  MAR 
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VII.    CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FURTHER  STUDY 

The  design  and  execution  of  this  experiment  has  resulted  in  several 
noteworthy  findings  and  recommendations  for  workers  pursuing  further 
development  in  this  area  of  numerically  integrating  conventional  and 
satellite  data . 

Based  on  the  availability  of  satellite  data  and  the  successful  pro- 
cedures developed,  (as  outlined  in  section  IV),  it  is  concluded  that  satel- 
lite data  inputs  to  numerical  analysis  is  operationally  feasible  and  worth- 
while.    Furthermore,  the  modifications  should  be  done  by  someone  (or  a 
group)  who  is  familiar  with  the  structure  and  behavior  of  both  the  numeri- 
cal analysis  scheme  and  the  PE  model.    To  retain  the  proper  time  and 
spatial  continuity  of  developing  systems  ,  the  modification  procedures 
should  be  pursued  on  a  continuous  basis,  not  just  "now  and  then." 
Likewise,  the  technique  should  be  applied  over  the  entire  numerical  grid 
and  not  just  in  one  area. 

The  completely  objective  SIRS  values  should  be  either  automatically 
processed  along  with  conventional  data  or  be  assigned  the  first  priority 
for  bogus  input  and  be  weighted  nearly  the  same  as  conventional  data. 
The  quasi-objective  ATS  winds  near  300  mb  and  the  Nagle-Hayden  tech- 
nique    500  mb  height  values  for  spiral  cloud  centers  are  reliable  for  use 
in  the  sparse  data  regions  and  should  be  treated  with  second  priority  for 
bogus  input o    The  other  subjective  pattern-recognition  techniques  can 
produce  a  large  quantity  of  approximate  500  mb  height  values  that  are 
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useful  for  sparse  data  regions  and  should  be  considered  as  tertiary 
priority  for  bogus  input. 

It  is  recognized  that  the  ultimate  effect  of  the  satellite-improved 
MOD  ANALS  to  generate  conclusively  favorable  prognoses  at  the  surface 
and  500  mb  was  not  achieved.     However,  the  MOD  PROGS  cannot  be 
considered  conclusively  unfavorable  either  in  view  of  the  fact  that  there 
is  definite  time  and  space  continuity  in  the  prognosis  of  those  pressure 
systems  introduced  in  the  analysis  through  satellite  observations, 
although  such  systems  were  significantly  damped  early  in  the  prognostic 
process.     Further,  the  merits  of  the  satellite-generated  modifications 
may  be  heavily  screened  by  behavioral  characteristics  of  both  the  initial- 
ization procedures  and  the  PE  model. 

It  is  also  noteworthy  that  the  incorporation  of  a  considerable  amount 
of  carefully  estimated  bogus  values  did  not  result  in  the  FNWC  PE  prog- 
nostic model  generating  any  spurious  errors  a  considerable  distance  down- 
stream of  the  modified  analysis  area. 

Finally,  the  reader  is  cautioned  against  making  sweeping  general- 
izations from  the  limited  results  reported  on  here.     Only  further  testing 
will  determine  the  universality  of  conclusions  as  drawn  from  the  small- 
sample  case  study  presented. 
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APPENDIX  A 
COMPUTER  RECTIFIED  PHOTOGRAPHS  OF  ATS-1  OBSERVATIONS; 
1800  GMT  9-12  MARCH  1971 
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APPENDIX   B 
COMPARISON  OF  NORTHERN  HEMISPHERE  SIRS    500  mb  HEIGHTS 
AND  FNWC  MONTEREY  500  mb  ISOHYPTIC  ANALYSES 
FOR  9-14  MARCH  1971 

During  the  period  9-14  March  1971,  SIRS  500  rnb  height  data  cover- 
age over  the  North  Pacific  Ocean  was  insufficient  to  give  a  perspective  of 
the  accuracy  and  potential  value  of  this  data  source.    To  provide  a  better 
study,  all  the  Northern  Hemisphere  SIRS  500  mb  values  for  9-14  March 
were  plotted  on  blank  charts,  overlain  on  the  corresponding  FNWC  500  mb 
operational  contour  analysis  maps,  and  the  FNWC  contour  height  compared 
with  the  SIRS  height.    This  approach  differs  considerably  from  the  study  of 
Nilsestuen     [161     where  the  SIRS  soundinjs  were  compared  with  actual  radio- 
sonde profiles  from  the  same  time  and  vicinity.     Since  a  SIRS  value  does 
represent  a  sampling  of  a  large  area  (a  square  120  NM  on  a  side),  a  com- 
parison with  the  500  mb  contour  analysis  is  perhaps  more  valid  than  that 
with  a  radiosonde  height  which  is  a  discrete  air  sample,  microscopic  in 
size. 

SIRS  soundings  north  of  70N  are  classified  by  NESS  as  "unreliable," 
and  used  as  an  input  at  NMC  with  lesser  "weight"   (analysis  significance). 
SIRS  soundings  south  of  7 ON  are  given  the  same  reliance  as  conventional 
soundings  in  the  NMC  objective  analysis.    This  study  considered  SIRS 
values  north  of  70N  separately  to  give  some  measure  of  their  reliability. 
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Two  tolerance  checks  were  tabulated  in  this  study:     1)  Category  1  — 
greater  than  30  m  difference,  and   2)  Category  2  —  greater  than  60  m  differ- 
ence.    SIRS  sounding  values  which  fell  into  the  second  category  (>   60  m) 
are  included  in  the  first  category  (  >  30  m). 

Table  III  sums  up  the  accuracy  of  SIRS  soundings  when  compared  with 
the  conventional  data  analysis,  on  a  hemispheric  basis.     For  the  eleven 
map  times  (09/0000Z  -  14/OOO0Z  March  1971),  an  average  of  56  soundings 
per  synoptic  time  went  into  the      update        NMC  500  mb  analysis.    The 
average  number  of  conventional  soundings  at  synoptic  map  time  is  approxi- 
mately 5  75;  so,  the  NMC    use    of  SIRS  data  represents  an  increase  of  10 
per  cent  coverage .    As  to  its  accuracy,  this  brief  study  showed  that  an  aver- 
age of  five  of  the  56  reports  differed  from  the  FNWC  conventional  data 
analysis  values  by  more  than  60  m.    The  modification  effort  made  by  the 
authors  in  the  North  Pacific  supports  a  tentative  conclusion  that  the  accu- 
racy of  these  differing  reports  is  in  favor  of  the  SIRS  data. 

The  results  for  the  polar  region  (north  of  70N)  lend  minor  credence  to  the 
present  NESC  classification  of  "unreliable"  for  the  SIRS  data.    An  average 
of  two  out  of  eleven  reports  were  greater  than  60  m  different  from  the  FNWC 
analysis  values.    However,  considering  the  scarcity  of  reliable  conventional 
reporting  stations  in  the  polar  region,   11  SIRS  soundings  represent  a  major 
increase  in  coverage. 
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APPENDIX   C 
COMMENTS  ON  PROCEDURES  USED  TO  DETERMINE  BOGUS  DATA 

Pursuant  to  the  desired  numerical  MOD  ANAL  is  the  development  of 
procedures  for  utilizing  all  available  options  for  obtaining  a  numerical 
MOD  ANAL  which  is  a  duplicate  of  the  manually-derived  MOD  ANAL. 
Bogus  values  of  pressure  or  height,  and  wind  to  modify  existing  centers, 
troughs  and  ridges, can  be  positioned  at  any  geographic  coordinates  and 
injected  into  the  analysis  scheme  at  sea  level  (100  points)  and  500  mb 
(50  points),  with  priority  over  most  rejection  criteria.     (A  greater  number 
of  bogus  values  will  be  acceptable  at  other  levels  in  the  FNWC  program 
in  the  near  future.) 

Bogus  values  that  blend  without  conflicting  with  the  conventional 
data  do  produce  good  numerical  MOD  ANALS.     If  the  desired  MOD  ANAL 
values  were  not  achieved  on  the  first  attempt,  the  pressure  or  height 
values  assigned  to  the  bogus  positions  in  the  unacceptable  area  required 
a  compensation  called  "overbogus"  which  is,  in  effect,  applying  the 
necessary  adjustment  in  pressure  or  height  to  force  the  analysis  toward 
the  desired  MOD  ANAL  value  as  it  smooths  the  data. 

When  the  "overbogus"  technique  was  not  successful,  it  was  neces- 
sary to  "override"  the  erroneous  conventional  data  by  saturating  a  specific 
area  with  numerous  bogus  values,  which  was  uneconomical  in  view  of  the 
limit  set  on  the  number  of  bogus  data  points  allowed  for  this  study. 
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Particular  cases  were  encountered  where  neither  one  nor  both  of  these 
techniques  used  together  was    sufficient  to  achieve  the  desired    MOD 
ANAL  values  over  the  area  of  concern  due  to  the  limit  of  points ,  inex- 
perience of  the  authors   and   design  criteria. 
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